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Abstract 

In this paper we study the locomotion of a shape-changing body swimming 
in a two-dimensional perfect fluid of infinite extent. The shape-changes are 
prescribed as functions of time satisfying constraints ensuring that they result 

1 I from the work of internal forces only: conditions necessary for the locomotion 

to be termed self-propelled. The net rigid motion of the body results from the 

O exchange of momentum between these shape-changes and the surrounding 
fluid. 

The aim of this paper is three-folds: First, it describes a rigorous frame- 
work for the study of animal locomotion in fluid. Our model differs from 

a previous ones mostly in that the number of degrees of freedom related to 
, , the shape-changes is infinite. The Euler-Lagrange equations are obtained by 

applying the least action principle to the system body-fluid. The formalism 
of Analytic Mechanics provides a simple way to handle the strong coupling 
^ between the internal dynamic of the body causing the shape-changes and the 

dynamic of the fluid. The Euler-Lagrange equations take the form of a cou- 
pled system of ordinary differential equations (ODEs) and partial differential 
00 equations (PDEs). The existence and uniqueness of solutions for this system 

are rigorously proved. 

Second, we are interested in making clear the connection between shape- 
changes and internal forces. Although classical, it can be quite surprising to 
select the shape-changes to play the role of control because the internal forces 
^■—^ they arc due to seem to be a more natural and realistic choice. We prove that, 

^ when the number of degrees of freedom relating to the shape-changes is finite, 

• ^ both choices are actually equivalent in the sense that there is a one-to-one 

relation between shape-changes and internal forces. 
^ Third, we show how the control problem, consisting in associating with 

each .shape-change the resulting trajectory of the swimming body, can be 
analysed within the framework of geometric control theory. This allows us 
to take advantage of the powerful tools of differential geometry, such as the 
notion of Lie brackets or the Orbit Theorem and to obtain the first theoretical 
result (to our knowledge) of control for a swimming body in an ideal fluid. 
We derive some interesting and surprising tracking properties: For instance, 
for any given shape-changes producing a net displacement in the fluid (say, 
moving forward), we prove that there exists other shape-changes arbitrarily 
close to the previous ones, that leads to a completely different motion (for 
instance, moving backward): This phenomenon will be called Moonwalking. 
Most of our results are illustrated by numerical examples. 
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1 Introduction 

1.1 History 

In the last decade, much work have been done by mathematicians to better under- 
stand the dynamics of swimming in a fluid. This interest has grown from the obser- 
vation that fishes and aquatic mammals evolved swimming capabilities far superior 
to those achieved by human technology and consequently provide an attractive 
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model for the design of biomimetic robots. Such swimming devices propelled and 
steered by shape-changes would be more efficient, stealthier and more maneuverable 
than if propeller-driven. This explains why, for instance, autonomous underwater 
vehicles are catching the attention of the petroleum industry for their possible use 
in the maintenance of off-shore installations. In the field of nano-technology, the 
design of nano-robots able to perform basic tasks of medicine is a challenge for the 
forthcoming years. 

Significant contributions to the understanding of the biomechanics of swimming 
have been made by Lighthill [!!)], Taylor [■->'), 'M>] Childress [ ] and Wu [38]. An 
interesting survey on the general theme of fish locomotion written by Sparenberg 
[33] is worth being mentioned as well. 

Experiments have shown that the vortices generated by the tail fins of fish play 
a crucial role in their locomotion and some models incorporate artificially produced 
vortices [23, 39, 37]. If we do not neglect the viscosity effects, the relevant model 
incorporates the non-stationary Navier-Stokes equations for the fluid coupled with 
Newton's laws for the fish-like swimming object. This perspective is adopted by 
Carling, Williams and Bowtell in [3], Liu and Kawachi in [20], Galdi in [s] or San 
Martin, Scheid, Takahashi and Tucsnak in [3(J]. However, and contrary to some 
common beliefs, forces and momenta acting on the fish body by shed vortices are 
not solely responsible for the net locomotion and among numerous mathematical 
articles studying fish locomotion, most of them address the case of a potential flow 
which is by definition vortex- free: let us mention here the works of Kelly and Murray 
[I-")], Kozlov and Onishchenko [!(>], Kanso, Marsden, Rowley and Melli-Huber [14], 
Melli, Rowley and Rufat [25] and Munnier [27, 28]. This is also the point of view 
we have chosen in this article. 

Although crucial for the design of autonomous underwater vehicles, results on 
control or on motion-planning for this kind of problem are very few, most of them 
focus on articulated bodies as in the works of Alouges, DeSimone and Lefebvre [2] 
(dealing with a three spheres mechanism swimming in a viscous fluid) or those of 
Mason [22] and Melli, Rowley and Rufat [24]. More authors have considered the 
problem of controlling immersed rigid solids. Let us mention for instance the paper 
[31] of San Martin, Takahashi, and Tucsnak, in which the control is chosen to be 
the relative fluid's velocity (thrust) on the solid's boundary, while in the paper [ ] of 
Chambrion and Sigalotti, the control is the impulse of the fluid also on the body's 
boundary. Finally, in [Ki, 17], the authors examine how a body with a rigid hull 
wich can modify the balance of its internal mass has the ability to steer and propel 
itself in a perfect fluid. 

The shape-changing body we consider in this paper (sometimes called amoeba 
for its similarity with this single-celled animal) is inspired by that of Shapere and 
Wilczek introduced in [.')2] and further discussed in [G]. The Euler-Lagrange equa- 
tions of motion are obtained following the method described in [l-l] for articulated 
rigid bodies (i.e. shape-changing bodies made of rigid solids linked by hinges) or in 
[27, 28] for more general deformations but adapted here to our inflnite dimensional 
model. Notice, however, that the main idea, consisting in invoking the least action 
principle of Lagrangian Mechanics for the overall system fluid-body, goes back to 
the works of Thomson, Tait and Kirchhoff in their studies of the motion of rigid 
solids in a perfect fluid. In his book [18, chap. VI, page 160], Lamb explains that: 
The cardinal feature of the methods followed by these writers consists in this, that 
the solids and the fluid are treated as forming together one dynamical system, and 
thus the troublesome calculation of the effect of the fluid pressures on the surfaces 
of the solids is avoided. 

The shape-changes are prescribed as functions of time and used as controls to 
propel and steer the amoeba. Since we want the motion to be self-propelled, these 
deformations have to result from the work of internal forces and torques only. It 
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entails some physical constraints previously discussed in [30, 27, 28]. Obviously, the 
fluid domain changes along with the shape of the amoeba. As in [23] and mostly 
in [.Ti], these changes are described by means of conformal mappings allowing the 
explicit computation of the fluid potential. Indeed, one of the main difficulties we 
are faced with in studying flsh-like locomotion is the precise analysis of the fluid 
potential with respect to the variations of the fluid domain. 

1.2 Main results 

The main results of this paper adress, on the one hand, the modelling and the well- 
posedness of the Euler-Lagrange equations and, on the other, the associated control 
problem. 

Our model for a swimming shape-changing body, thoroughly described in Sec- 
tion 2, constitutes the first novel concept set out in this paper. Although profoundly 
inspired by [6], it has been substantially improved. In particular, in the article of 
Shapere and Wilczek, the inertia of the amoeba is neglected whereas here, the 
time-evolving mass-distribution inside the animal is taken into account. Further, 
the shape-changes we consider are richer than those of the model in [32] , for they 
have an infinite number of degrees of freedom. This latter improvement leads us, 
in Proposition 3.1, to extend the validity of Kirchhoff's law to the case where the 
decomposition of the fluid potential into a linear combination of elementary poten- 
tials involves an infinite number of terms. A generalization of the classical notion 
of mass matrix is also required and this task is carried out in Subsection 3.2. The 
elementary potentials are solutions of boundary value problems set on the fluid 
domain which changes along with the shape of the amoeba. We prove in Theo- 
rem 3.1 that the elementary potentials, seen as functions of the shape of the body 
and valued into suitable Sobolev spaces, are smooth. As a straightforward con- 
sequence, we deduce in Theorem 3.2 that the mass matrices and the Lagrangian 
function, seen also as functions of the shape of the amoeba and valued into suitable 
spaces of bilinear mappings, are also smooth. This result is used in Subsection 3.4 
to derive the Euler-Lagrange equation (3.33), a system of ODEs that governs the 
dynamic of the system fluid-body. It is given in a form very convenient to study 
locomotion problems for it gives the rigid motion of the body with respect to the 
shape-changes. The sharp regularity results obtained for the mass matrices are also 
required to prove the well-posedness of this system of ODEs, see in Proposition 3.3. 

We next study the connection between shape-changes and internal forces. In 
Section 3.7, we show how the internal forces causing a given shape-change can be 
computed. In Theorem 3.3, we prove that if we consider only finite dimensional 
deformations, there is a one-to-one relation between shape-changes and internal 
forces: in other words, it is equivalent, in this case, to select for controlling the 
swimming animal either the shape-changes or the internal forces. 

The usual control problems associated with our dynamical system could be 
stated as follows: 

• Controllability: Is it possible to find shape-changes that propel the swim- 
ming animal from a given starting position to a specific end one? 

• Tracking: Is it possible to find shape-changes that allow the amoeba to follow 
(approximately) any given trajectory? 

We show that the answer to these two questions is positive and actually we prove 
more. In Theorem 4.1, we claim that the swimming animal can not only follow 
approximately any prescribed trajectory, but also while undergoing approximately 
any prescribed shape-changes. This surprising result is, to our knowledge, the first 
theoretical controllability result obtained for a body swimming in a perfect fiuid. 
Most of these results are illustrated by numerical examples in Section 5. 
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1.3 Outline of the paper 

The modeling is performed in Section 2, which deals mainly with the description of 
the shape-changes and the kinematics of the problem. Dynamics is treated in Sec- 
tion 3 where we derive the Euler-Lagrange equation. In Subsection 3.7 we discuss 
the equivalence between controlling by shape-changes and by internal forces. Sec- 
tion 4 is devoted to the control problem: the main controllability results are stated 
(with comments) in Subsection 4.1 and then proved in the remaining Subsections. 
Finally, Section 5 is dedicated to the study of trajectory-planning and to numerical 
simulations. 

Appendices A provides additional material used in Sections 2 and 3 and Ap- 
pendix B contains classical results of geometric control theory used in Section 4. 

2 Setting of the problem 
2.1 Notation 

In this Subsection, we introduce the main notation and the functional framework. 

2.1.1 Systems of coordinates 

Let (ei, denote a reference Galilean frame by which we identify the physical space 
to R^. At any time the amoeba occupies an open smooth connected domain A and 
we denote by :— \^ the open connected domain of the surrounding fluid. The 
coordinates in (ei,e2) are denoted with lowercase letters x = {xi,X2)'^ (the super- 
script T standing for the matrices or vectors transpose) and are commonly called 
the spatial coordinates (see for instance [21, chap 15]). For any x = (xi, X2)'^ G R^, 
a;"*" := {—X2,xi)'^ stands for the vector x positively quarter turned. 

Attached to the microorganism, we define also a moving frame (e|,e2). We 
choose it such that its origin coincides at any time with the center of mass of 
the swimming animal and we denote by x* = (x^jXj)^ the related so-called body 
coordinates. In this frame and at any time the amoeba occupies the region A* and 
the fluid the domain J^* R^ \ A* (see flgure 1). More generally, quantities with 
an asterisk are expressed in the moving frame. 

We define also the computational space endowed with the frame (Ei,E2) and 
in which the coordinates are denoted z — (zi, 22)^. In this space, D is the unitary 
disk and := R2 \ D. 

Throughout this paper, we will use the same notation n to represent the unitary 
normal to dA — dT or dA* — dT* directed toward the interior of the amoeba. 
We will sometimes use complex analysis and identify the space we are working in 
with the complex field C. In this case the notations introduced above will turn 
into X = Xi + 1x2, X* — x\ + 1X2 or z = Zi + iz2 with = — 1 and x^, x^,, zj, £ R 
(k = 1,2). The complex conjugate of z = zi +iz2 is z = zi — iz2. We will sometimes 
also mix this notation with the polar coordinates (r, 9) £ R-(- x R/27r, z — re'^. 

2.1.2 Function spaces 

Let E and F be Banach spaces and assume that X is a compact subset of E. The 
vector space C'^{K,F) (m an integer, m > 1) of the functions m times continu- 
ously differentiable from K into is a Banach space once endowed with the norm 
of W"'°°(if, F) (uniform convergence in K of the function and all of its partial 
derivatives up to the order m). 

Weighted Sobolev and Lebesgue spaces required in the resolution of boundary 
value problems are introduced in the Appendix, Subsection A. 2. 
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2.1.3 Multilinear, polynomial and analytic functions 

Let El, . . . ,Ek {k > 1) be Banach spaces. The set consisting of all the continuous, 
fc-linear mappings from E Ei x . . . x Ek into F is denoted Ck{E,F) (we will 
drop the subscript k when fc = 1). It is a Banach space whose norm is classically 
defined by: 

l|A||£,(ij,F) := sup |l(A,ei,...,efe)|i;^, {KeCk{E,F)). 

(ei,.-.,efc)eE, 

ll=jllEj.=l. (j = l, •••,'!) 

When = R, we are dealing with multilinear continuous forms and we denote 
merely Ck{E) := Ck{E,F). 

We call polynomial function from a Banach space E (or from only a subset of 
this space) into a Banach space F any function P such that there exists an integer 
p G N (the degree of the polynomial), Aq S F and p mappings Ak g Ck{E^,F) 
(k — 1, . . . ,p) such that: 

p 

P(e) := Ao + ^(Afc,e,...,e), Ve e 

fe=i 

We denote 'P{E, F) the set of all the polynomial functions from E into F. Observe 
that in particular C{E,F) C V{E,F). We easily prove that if Ei, E2 and E^ are 
three Banach spaces and Pi € P(i;i, £'2), P2 e V{E2, E3) then P20-P1 e 'P(£'i, £^3). 

We call analytic function from a Banach space E into a Banach space F any 
function / such that there exists R> (the radius of convergence), Aq ^ F and a 
sequence (v4fc)fe>i with A^ € Ck{E'',F) satisfying: 

J2 \M'\\Mc,(EKF) < VA e R, |A| < i?, 

k>l 

and 

/(e) := Ao + ^(Afc,e,...,e), Vee£, Me < R.. 

k>l 

We refer to the book [ , §4], for further details on analytic functions in Banach 
spaces. 



2.1.4 Banach spaces of series 

We denote any complex series by c := (cfc)fe>i where for any fc > 1, Cfc := Ofe + ibk, 
ak,bk S R. Most of the complex series we will consider in this article live in the 
Banach space: 

S := {{ck)k>i e : ^ fcdflfcl + \bk\) < +00}, 

k>l 

endowed with its natural norm ||c||5 := X]fc>i ^(l^^fcl + l^fcD- The unitary ball of S 
is denoted B and we will require the following open subset of S, containing B: 

V := |c := (cfc)fe>i e 5 : sup I V(fc + l)ck+iz^ < l] 
^ zedD I ^ > 

Let us introduce also the Hilbert space: 



fc>0 



k>l 

whose natural norm is denoted ||c||r- From the obvious identity ||c||r < ||c||s, we 
deduce that S cT. Further elementary results about these spaces are given in the 
Appendix, Subsection A.l. 
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2.2 Kinematics of the shape-changing amoeba 



2.2.1 Description of the shape-changes 

The shape-changes of the amoeba are described with respect to the moving frame 
(01,62) by a diffeormorphism x(c), depending on a shape (or control) variable 
c E T), which maps the unitary closed disk D of the computational space onto the 
domain A* of the physical space. We can write, according to our notation, that for 
any c E T), A* = x(c)(i5) and x* = x(c)(z), {z E D). The map x(c) is defined for 




Figure 1: Quantities are denoted with an asterisk when they are expressed in the 
moving frame (e*). 



all c e P by: 

X(c)(z) :=xo(z) + (xi,c)(z), (zed), (2.1) 

where Xoi^) ^ and (xi,c)(z) := J2k>i'^kz'^' ^ Introducing in polar 

coordinates {r,0) E [0, 1] x R/27r the vectors fields: 



cos{k9) 
-sm{k9) 



and C/^(r,( 



'sm{kd) 
cos{k9) 



(2.2a) 



we have the equivalent definition: 



X{c){r,6) r 



cos{9) 
sin{9) 



J2akU^{r,0) + bkUj:ir,9). (2.2b) 



k>l 



We next introduce likewise the function (f){c) that maps onto J-*. It is defined 
for all c e I? by: 

0(c)(z) :=z + ^CfcZ-^ (zEn). (2.3) 



fc>i 



Since z = 1/z for all z E dD, we deduce that x{^)\dD = 'P{c)\dn and the following 
mapping is continuous in C for all c e P: 



X{c){z) a ZED, 
0(c)(z) if z eCI. 



^c)iz) := 

We can summarize the main properties of x and (f> in the following proposition: 
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Proposition 2.1. For all c eV, x(c) : D ~^ A* and 0(c) : Cl T* are both well- 
defined (the series in (2.1) and (2.3) converge for all z in D and Cl respectively) and 
invertible. Further, x(c)|d is a diffeomorphism, (/)(c)|q is a conformal mapping 
and <i>(c) is an homeomorphism from C onto C. 

Proof. Let us denote I?x(c)(z) the Jacobean matrix at a point z G I? of x(c)(z), seen 
as a function from into R^. Identifying the matrix-vector product Dx{c){z)h 
for all h := {hi, /i2)"^ S R^ with its complex expression in which h — hi + i/i2, we 
can write that: 

Dxic){z)h=h+ (^(fc + l)cfe+iz'=)/i. (2.4) 

*;>0 

The series in the right hand side term is the conjugate of the holomorphic function 
z G D 1-^ X]/c>o(^ whose maximum is achieved, according to the maxi- 

mum principle, on dD. We deduce that for all c G T), x(c) — Id is a strict contraction 
in D. Invoking the local inversion Theorem, we deduce that for all c G I?, x(c) is 
a local diffeomorphism in D. 

We use roughly the same arguments to prove that (/)(c) is locally a conformal 
mapping. Starting from the expression (j)'{c){z) — l = (— 1/z^) J2k>oi^^^)^^k+iz^'^ , 
we deduce, according to the maximum principle and since z^'^ — z^ for all z £ dD, 
that |(/.'(c)(z)-l| < sup^ga^ |X;fc>i(fc + l)cfe+i^''| = sup,^9o\T,k>iik+l)ck+iz''\. 
It entails that for all c G V, 0(c) — Id is a strict contraction in fi and then that (j){c) 
is locally a conformal mapping in f2. 

Since both mappings x(c) — Id and (j){c) — Id are strict contractions, we draw 
the same conclusion for <&(c) — Id in the whole complex plane C. For any z £ C, 
the function z E C i-^ z — ($(c)(z) — z) being also a strict contraction, Banach 
fixed point Theorem ensures the existence and uniqueness of z g C such that 
z — (<i>(c)(z) — z) — z i.e., (f>(c)(z) = z. The proof is then completed. □ 



2.2.2 Body's volume, density, mass and inertia momentum 

From the relation x* :— x(c)(z), {z E D), we deduce that the area elements dx* 
and dz of respectively A* and D can be deduced one from the other by the identity: 

dx* := I det Dx{c){z)\dz, {z € D, x* := x{c){z)). 

This entails that the density p* of the deformed amoeba A* can be deduced from a 
given constant density po > by the conservation-of-mass principle: 

- Id.. z,,(c)(:(c)-.(.-))r '-''-■'''^ 

Lemma 2.1. For all c E V, the volume of the amoeba is: 

Yol{A) = 7r(l - ||c||2-). (2.6) 
Proof. Starting from the expression (2.4) of Dx{c), we obtain that: 

detDx{c){z) = 1 - (^fccfcz'=-i)(^A:cfeZ^-i), (z E D), 

k>l k>l 

and it is clear that the right hand side term is positive if c S P. Expanding it, we 
get in polar coordinates {r,9) E [0, 1] x R/27r: 



I det Dxic)ir, d)\ ^1-^.^2 kk' r''+''' -^a^ak' + hkhy) cos((fc - k')9) 

k>l k'>l 

+ kk'r''+'''-^{akbk' - ak'bk) sin((fc - k')9). (2.7) 



fc>l k'>l 
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Since the volume of the ameoba is given by Vol(.4) — | det Dx{c){z)\dz, straight- 
forward and easy computations lead to formula (2.6). □ 

We can now define the element of mass in D by dmo := padz, and likewise 
dm* := p^dx* , is the element of mass in A* . 

As the diffeomorphism x(c) is modeling physical shape-changes, it has to satisfy 
some constraints. Let us consider a continuous and piecewise control function 
t e [0,T] 1-^ c{t) {T > 0). This regularity entails that its time derivative c{t) 
exists for all but a finite number of t in [0,T]. We denote c = {ck)k>i ^ S and for 
all fc > 1, Ck := ctfe + ibk- We introduce as well: 

x(c)(z):^^(z) = (Xi,c), {zeD,t>0). (2.8) 

Because of the incompressibility of the fiuid, its volume has to be constant. We draw 
the same conclusion for the volume of the amoeba because its volume is nothing 
but the complementary of the volume of the fiuid. According to (2.6), it means that 
the function t e [0,r] i-^ ||c(t)||r has to be constant for the control function to be 
physically allowable. We introduce the notation /i := ||c(0)||r and the subset of S: 

£{p) := {ce5 : ||c||r = m}- (2.9) 

Remember that for the map x{c) to be injective, c has to belong to V. We define 
then also: 

£'{p) ■.= {ceV : ||c||r-/i}. (2.10) 
According to Parseval's identity: 

E k\ck\' < E ^^'l^^l' =i r ( + l)ck+ie-^''yde 



k>l k>l k>Q 



< sup I V(A:-H l)cfe+iz'= 



2 



and we deduce that £*{p) is non empty if and only li p < 1. The constant volume 
of the amoeba can next be rewritten: 

Vol(^) = 7r(l - /x^), {t>Q). (2.11) 

Differentiating with respect to time identity (2.6), we get an equivalent formulation 
for the conservation of the amoeba's volume: 

^k{akak + hkbk) = Q, (^ > 0), (2.12a) 

fe>i 

or, with notation of Subsection A.l of the Appendix: 

(G(c),c)=0, (t>0). (2.12b) 

Remark 2.1. Quite surprisingly, the constraint on the body's volume is not required 
any longer in dimension 3 as proved in [ ]. A physical explanation is that the 2d 
case can be seen as a 3d model in which we consider the swimming animal as an 
infinite cylinder of section A. Hence, any shape-change of the body which does not 
preserve its volume, entails an infinite variation of the fluid's volume - which is 
impossible because the fluid is incompressible. With a real 3d model, although the 
fluid may be also incompressible, the finite variations of the fluid's volume due to 
the shape-changes can, in some sense, be neglected when compared with the infinite 
overall amount of fluid. From a mathematical point of view, the boundary value 
problem driving the motion of the fluid requires condition 2.12a (or 2.12b) to be 
well-posed in dimension 2 and not any longer in dimension 3. 
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The mass of the amoeba is: 

m := dm* = / dmg — irpQ. (2-13) 
J A* J D 

If we assume that the body is neutrally buoyant, we get the equality p/Vol(^) = m 
where pf > is the given constant density of the fluid and we deduce that the 
densities pf and po are linked by the relation: 

Pf{l~p')=p^. (2.14) 

The inertia momentum depends on the shape of the amoeba and reads /(c) :— 
Xa- \x*\^ dm* , or equivalently upon a change of variables /(c) :— |x(c)(z)P dmo. 
It can also be easily computed in terms of the control variable: 

/(c) = ^Po(^ + EfcXTl^^-l')- (^-^^^ 

fc>l 

In this form, we sec that /, seen as a function of c valued into R, belongs to V{'D, R). 
2.2.3 Constraints where the motion is self-propelled 

Since we assume that the shape-changes of the amoeba are produced by internal 
forces and torques only (by definition of a self-propelled motion), Newton's laws 
ensure that the linear and the angular momenta of the animal with respect to 
its attached frame (e|,e2) have to remain unchanged when it undergoes shape- 
changes. Hence, we get the condition: 



d 
di 



^ j x*dm*^ =0, {t> 0), 
which can be rewritten, upon a change of variables and taking into account (2.5): 

x(c)dmo = 0, (t>0). (2.16) 
For the angular momentum, the same arguments yield: 

x(c)-x(c)^dmo = 0, (t>0). (2.17) 



D 

Condition (2.16) is actually intrinsically satisfied for any control function t ^ c{t). 
For condition (2.17), observe first that, with definition (2.2), {Uj^)-^ = C/^ and 
(U^)'^ = —U^, which leads after some basic algebra to the identity: 

J2]A^(^kak-dkbk) = 0, {t>0), (2.18a) 

k>l 

or equivalently, with the notation of Subsection A.l: 

(F(c),c)=0, {t>0). (2.18b) 

This constraint together with (2.12b) lead us to introduce the notion of allowable 
control: 

Definition 1 (Physically allowable control). A smooth function t e [0,T] i-^ c{t) e 
S (for some real positive T) is said to be Physically allowable when: 

• There exists p > such that c{t) G £'{p) for all t > 0. 
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• Constraint (2.18b) is satisfied for all t £]0,T[. 





ij 


^1 





Figure 2: Examples of physically allowable shape-changes (/j, = 1/^/12, \\c\\s < 1). 
The colors gives the value of the density inside the animal (neutrally-buoyant case). 



2.3 Rigid Motion, Velocity 

The overall motion of the amoeba in the fluid consists in the superimposition of 
its shape-changes with a rigid motion. The (prescribed) shape-changes have been 
described in the preceding Subsection along with the (unknown) net rigid motion 
results from the exchange of momentum between the shape-changes with the sur- 
rounding fluid. It is described by elements q := {r,9)^ of Q := x R/27r where 
r :— (ri,r2)'^ € is a vector giving the position of the center-of-mass of the body 
and 9 € R/27r an angle giving its orientation with respect to (61,62); see Fig. 1. 
If we denote by R{0) G S0(2) the rotation matrix of angle 9, then we have the 
relations R{9)ej = e* (j = 1,2). 

Let the shape-changes be frozen for a while and consider a physical point at- 
tached to the body undergoing only a rigid motion. Then, there exists a smooth 
function t G [0,T] ^ q(t) ~ {r{t),9{t)) G Q (T > 0) such that the point's co- 
ordinates in (61,62) be given by x = R{9)xo + r (a;o € R^ being the coordinates 
at the time t = 0). Next, compute the time derivative of this expression and 
denote by q :— {r,u!) G R^ the time derivative of q. Since we classically have 
dgR{9)R{9)^x — for all x — (xi, ^2)'^ G R^, we deduce that the Eulerian veloc- 
ity of the point is Vr{x) = uj{x — r)-^ + r. It can also be expressed in the moving 
frame (6^,62) and reads v*(a;*) = uj^x*)-^ + r* where r* := R{9)'^v. This leads us 
to introduce also the notation q* := (r*,a;)-'" G R'^. 

Let us return to the general case where the shape-changes are taken into account. 
The coordinates in (61,62) of a physical point attached to the amoeba are given by: 
x = R{9)x{c){zq)+v, where at the time t = 0, c(0) = Cq and x(0) — xo — x{^o){zo) 
for some zq G D. Observe that we can always assume that at the time t — Q, 
q = qo := (0, 0)-^. We deduce that the Eulerian velocity at a point a: of ^ is: 

v(x) = u;{x - r)^ + r + R{9)x{c)[x{c)-\R{9f {x - r)))]. 

When expressed in the moving frame (6^^,62) it reads: 

v*{x*)^{ujx*^+r*) + x{c){x{c)-Hx*)), {x*gA*)- (2.19) 

2.4 Potential flow 

The fluid is assumed to be incompressible and inviscid. We denote pf > its 
constant density and we set dm^ :— pjdx* the element of mass in J-'* and dm^ := 
Pfdz the element of mass in fi. We seek the Eulerian velocity u* of the fluid, 
expressed in (6^,62), as the gradient of a potential function ip: 

u*:^Vip inr*. (2.20) 
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The incompressibility of the fluid entails that V • u* = and hence: 

Aip*^0 inT*. (2.21a) 

The classical non penetrating (or slip) condition for inviscid fluid leads to the equal- 
ity: 

u* • n = V* • n on dT* , 
and yields the following Neumann boundary condition for ip: 

dn(p — V* ■ n on dJ-* . (2.21b) 

The boundary value problem (2.21) admits a weak (or variational) solution in the 
weighted Sobolev space H\i{T*) (see Appendix, Subsection A. 2). The potential 
function is actually only defined up to an additive constant. It does not matter 
since we are only interested in V(/5 which is uniquely determined. Note that the 
potential function does depends on both c (linearly through the boundary data) 
and c (through the domain T* ) . 

2.5 Lagrangian of the system fluid-amoeba 

Since we neglect gravity, the Lagrangian function reduces to the kinetic energy of 
the system fluid-body. Because of relations (2.16) and (2.17), there is a decoupling 
between the kinetic energy of the body due to its rigid motion and that due to its 
shape-changes: 

:= + \\^?m + |x(c)(x(c)-^(x*))pdm*, 

the last term being the kinetic energy of deformation. It can be further computed 
as follows: 

\x{c){x{c)-\x*))\^dm*^ [ |x(c)(z)pdmo = ^poEM-- (^.22) 

•^^ k>i 

On the other hand, the kinetic energy of the fluid reads: 

■= lul'dm} \V^\^dm). (2.23) 

The Lagrangian function of the system fluid-amoeba next reduces to: 

L:=K^ + K^, (2.24) 

and turns out to be a function of (q*, c, c) e (R^ x R) x I? x 5. More precisely, for 
any fixed c G "D, L{c) is a quadratic form in (q*,c). It is worth remarking that it 
does not depend on r and 6 due to the isotropy of our model with respect to the 
position and orientation of the body in the fluid. 

Remark 2.2. As already mentioned, if we do not neglect gravity but rather assume 
that the body is neutrally-buoyant we have to add to our model the relation (2.14) 
linking pf and po. However, (2.14) only ensures that the upthrust is null. The 
torque applied on the amoeba by the buoyant force is equal to: 

where g is the standard gravity and v f is the center of buoyancy given by: 

Vf := — / dmr. 
mJA 
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The torque is null only when r f — r what is not always verified under solely as- 
sumption (2.14). 

The main consequence of taking into account buoyancy is to lose the isotropy 
of our model (as seen by an observer attached to the body and, without buoyancy, 
all of the positions and directions in the fluid are equivalent). In this case, the 
Lagrangian function depends also on q and the Euler- Lagrange equations turn out 
to be much more involved. Models with buoyancy, several swimming bodies and 
bounded or partially bounded fluid domains are studied in [28]. 

3 Euler-Lagrange equations 

The aim of the Section is to obtain the Euler-Lagrange equation that governs the 
dynamics of our system. 

3.1 Elementary potentials 

KirchhofF's law states that the potential function can be decomposed into a linear 
combination of elementary potentials, each one being associated with a degree of 
freedom of the system (which are here: the translations of the body along (j = 
1,2), the rotation and all of the elementary shape-changes governed by the variables 
Cfc {k > 1)). This law is classical when the number of degrees of freedom is finite 
but it must be adapted to our infinite dimensional model. This is the first goal of 
this subsection. 

The elementary potentials, as being solutions of Neumann boundary value prob- 
lems set on the fluid domain, depend implicitly on the shape of the amoeba, i.e. 
on the control variable c. The second issue we will address in this subsection is 
to examine how smooth is this dependence. To carry out this task, we will use a 
conformal mapping (/)(c) that maps the fixed domain f2 (the exterior of the unitary 
disk) of the computational space onto the fluid domain of the physical space. 
Each elementary potential, once composed with (/)(c) will yield a function defined 
in the fixed domain f2 and whose dependence with respect to c is explicit. 

3.1.1 Definitions, Kirchhoff's law 

We begin with decomposing ip into the sum of ip^ (the potential associated with the 
rigid motion) and (p"^ (the potential associated with the shape-changes). Both are 
harmonic in J-'* and satisfy the following Neumann boundary conditions: dn<p^ = 
{ujx*-^ + r*) -n and dnp"^ = {x{c){x{c)~^ {x*))) ■ n on dJ^* . We wish now to obtain 
a further decomposition of ip"^ into an infinite linear combination (i.e. a series) of 
elementary potentials each one associated with an elementary shape-change. To 
this end, we introduce the functions (pf. and ip\ {k > 1), all of them being harmonic 
in J^* and which satisfy the following Neumann boundary conditions: 

dnipt:^U^{xicy\x*))-n-kak on dT* for all fc > 1, (3.1) 
dn^pl := Ul{x{c)-^{x*)) -n-khk on dT* for all k > 1. (3.2) 

The extra terms kak and kbk have to be added for the boundary data to satisfy 
a so-called compatibility condition necessary to ensure the well-posedness of the 
Neumann problems (see Subsection A. 2). From a more physical point of view, the 
elementary shape-changes driven by the shape variables ak and bk do not preserve 
the volume of the amoeba and has to be suitably modified. Observe however that 
under condition (2.12), we recover 

J2 akdn^t + hdnipi = a„/. (3.3) 

k>l 
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The question is now: does this identity entails the equahty X]fc>i o,kfk~^^k'Pk ~ f'^'^ 
When the number of terms in the sum is finite, it is nothing but classical KirchhofF's 
law. The question can hence be simplified into: what is the topology the left hand 
side series of functions in (3.3) has to converge for, that ensures the equality? And 
the answer is: the topology of L^(9jF*), because solutions of Neumann problems 
in H]^{!F*) depend linearly and continuously on their boundary data in L^{dJ-*), 
as detailed in the Appendix, Subsection A. 2. Actually, for any c £ T) and c £ S, 
we easily confirm that the series of functions in the left hand side of (3.3) converge 
normally on dT* and hence also in L^{dJ-*). 

To emphasize the dependence of ip'^ with respect to c and its linear dependence 
with respect to c, we denote it rather {ip'^{c),c). 

It remains for us to introduce a decomposition for the potential ip^ . So, let 
us define the elementary potentials (p'j (j = 1,2,3) as being harmonic functions 
in !F* satisfying the Neumann boundary conditions: dn^^ = ' (j = 1,2) and 
9ji(/?3 = (x*)-*- • n on dJ-* . We can now state: 

Proposition 3.1 (Potentials decomposition). For any allowable control (in the 
sense of Definition 1), we have the identities in H^{!F*): 

(p'' ^ rlifil + r^(p2 + i^iPs, (3.4a) 
(/(c),c)=^afe^^ + 6,/fc, (3.4b) 

k>l 

^ = ^'-+(/(c),c). (3.4c) 

Although it does not appear in the notation, the potentials ip^ {i — 1,2,3) and 
</j^, if \ {k > 1) do obviously depend on c since the domain JF* does. 



3.1.2 Dependence of the elementary potentials with respect to c 

We use complex analysis to compute the elementary potential functions so we iden- 
tify both the computational and the physical spaces with the complex plane C. 
As already mentioned, so as not to overload the notation we will mix the complex 
notation z = zi + iz2 or x* ~ x\ -\- («^ = —1) with the real one z = (zi, 22)"^, 
X* — {xl,X2)'^ and even with the polar coordinates {r,0), r — \z\ and 6 = Arg (z) 
(i.e. z = re*^) when necessary. Remember that D is the unitary disk of the com- 
putational space, U := C \ D and that for all c G D, the mapping 0(c) is defined 
by (2.3). We get the following expression for the unitary normal to dA*: 

n{x*) := n,{x*) + tn2{x*) = -z^^^^ {x* = 0(c)(z), z G dD), 

|(/)'(c)(z)| 

where 4>'{c) is the complex derivative of 0(c). We introduce next the functions $J(c) 
(j = 1,2,3) defined by: 

CJ(c)(z) := ^;:(x*), {x* = 0(c)(z), z e n). (3.5) 

Since (/)(c) is a conformal mapping, the functions CJ(c) are harmonic in fl and we 
compute that: 

^^^|^a„e;(c)(^) = a„^K^*)' {^*^^{c){z),zen). (3.6) 

The Neumann boundary conditions for ip"^ lead to the following boundary conditions 
for^; (j = 1,2,3): 

a„er(c)(^) = -5R(^0'(c)(z)), (3.7a) 
dnC2ic){z) = -^{z<f>'{c){z)), (3.7b) 



dnC3{c){z) = -3(0(c)(z)z0'(c)(z)), (z e dD). (3.7c) 
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We proceed likewise for the elementary potentials related to the shape-changes: We 



define ^^(c) and ^^.(c) for all A: > 1 by: 

Cdc){z) ■■= Vl{x*), (3.8a) 
ek{c){z):^v\{x*), (x* = 0(c)(z), z e f]). (3.8b) 

These functions are harmonic in U, and satisfy the Neumann boundary conditions: 

dn£.t{^){z) = -5R(z'=+V(c)(z)) - fcflfe, (3.9a) 
dnil{c){z) = -3(z'=+V(c)(z)) - kbk. [z e dD). (3.9b) 



It is clear, applying the results of Subsection A. 2, that all of these functions are 
well-defined in the weighted Sobolev space iJj^(ri). What we are interested in is to 
study their regularities, seen as functions of c valued in H^{fl). We invoke again 
the linear-continuous dependence of the solution in iJ^(n) of a Neumann boundary 
value problem with respect to its boundary data in L^{dD) and the problem is 
reduced to the study of the dependence of the boundary data (3.7) and (3.9) in 
L^{dD) with respect to c e P. Some simple estimated based on the identities: 

z'=+V'(c)(^) = z'+' - 5] (j- + kf-^, (3.10a) 

j>i-k ^ 

WMz<l>'{c){z) - 1 + ^ Cfcz'=+i - J2 kcjckz^-", (3.10b) 

fc>i fc>i fc>ij>i 

available for all z G dD, lead us to deduce that this dependence is polynomial. It 
entails: 

Lemma 3.1. The S^-type functions defined in (3.5) and (3.8), seen as functions of 
c valued in Hj^^fl) belong to V{'D, Hj^{fl)). 

However, we need to prove a little bit more. Indeed, let us define also by 
^'^{z) := (p'^{(t){c){z)) for all z G D. Again, to emphasize the dependence of with 
respect to c and c, we denote it rather (^'^(c),c). According to identity (3.4b), we 
get the decomposition in i7^(il), for all c S I? and all c G S: 

{^'ic),c) = J2 + ^'»^^('^)' (3.11a) 

k>l 

and we deduce that can be seen as a function of c valued in C{S, iJj^f (fi)). We 
prove in the Appendix, Subsection A. 3 that the function c e I? e -^(15, i?^(f2)) 
is also polynomial. 

Finally, let us introduce ^ and defined by ^(z) := Lp{<j){c){z)) and ^^{z) := 
Lp^{(l){c){z)). According to Proposition 3.1, we can state that, for all c e I? and all 
c e 5: 

^'' = m+r;^2+^C3, (3.11b) 

^ = r+(e'(c),c). (3.11c) 
We can now summarize all of the results obtained in this Subsection: 

Theorem 3.1. • Well-posedness: For any c E V and any c E S, the 

functions ^, (^'*(c),c), ^J(c) (j = 1,2,3^, ^^(c) and ^^(c) (k > 1) are well 
defined as elements of H^f^fl). 

• Decomposition: For any allowable control (c, c) eV x S, identities (3.11) 
hold. 
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. Regularity: G V{V,H],{n)) (j = 1,2,3;, G 'P{'D.H],{n)) (k > Ij 

and C'' e 

The properties of conformal mappings allow to write that 




So from now on we will refer to ^ as being the potential function in place of and 
likewise, we will call U = 1,2,3), ^g(c) and ^^(c) {k > 1), the elementary 

potentials. 



3.2 Mass matrices 

Mass matrix is a central notion in the modeling of fluid-structure interaction prob- 
lems. It can be deflned as the polarization of the kinetic energy of the system, seen 
as a quadratic form. Remember that in our case, the kinetic energy coincides with 
the Lagrangian function defined in Subsection 2.5 and that, for any fixed c g 2?, 
L(c) is a quadratic form in (q*,c) G R'^ x S. We define then M(c) as being the 
bilinear symmetric form on (R'^ x S) x (R^ x S) such that: 

L(c,q*,c) = i(M(c),(q*,c),(q*,c)). 

We next decompose it into M[''(c), a bilinear symmetric form on R^ x R^ (that can 
be identified with an actual 3x3 symmetric matrix) , N(c) a bilinear form on 5 x R^ 
and M'^(c) a bilinear symmetric form on S x S such that: 

(M(c), (c, q*), (c, q*)) = (M^'^(c), q*, q*) + {M''{c),c, c) + 2(N(c), c, q*). 

We are interested in finding explicit expressions for M'"(c), N(c) and M'*(c) and in 
studying their dependence with respect to the control variable c. 



3.2.1 Mass matrix related to the rigid motion 



We consider first 1^(0), the classical mass matrix of the amoeba associated to its 
rigid motion. The decomposition of the potential function given in (3.11b) leads us 
to introduce the symmetric 3x3 matrix: 



r(c) := 

m ■ 

771 

/(c) 



7^VC5:(c) • V^[(c) dmO . . . ^Ve[(c) . VC5(c) dmj 
.LVS(c) • Ver(c) dmj ■■■ /oVe5(c) . Ve5(c) dmO 



, (3.12) 



where we recall that m > is the constant mass of the amoeba and /(c) is its inertia 
momentum given in (2.15). The kinetic energy due to the rigid displacement of the 
amoeba can be written as the matrix-vectors product: (l/2)(r*,a;)M'"(c)(r*, w)"^. 
The latter matrix in the right hand side of (3.12) is usually referred to as an added 
mass matrix, relating here to the rigid motion of the animal. 

We deal now with the kinetic energy due to the shape-changes, considering 
separately the infinite and finite dimensional cases. 
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3.2.2 Infinite dimensional case 

To define the elements of the matrices N(c) and M''(c) we use the canonical basis 
{f ^, f ^, f'^} of R'^ and the Shauder basis {a-', h^,j > 1} of 5 defined in the Appendix, 
Subsection A.l. We have, for any 1 < fc < 3 and j > I: 

(N(c),a^f^•) := / V4^(c).V^;(c)dmO, (3.13a) 
Jn 

(N(c),b^f^•) / VC^(c) . Ve)(c)dm°. (3.13b) 
Denoting Sj. the Kronecker symbol, the entries of M''(c) are, for all j, k > 1: 

(M'^(c),a^a'^-) ^ VC;(c) • VC,"(c)dm« + (3.14a) 

(M^(c),b^b^■> := Ve)(c) • Ve^(c)dm° + (3.14b) 

(M'^(c),a^b^-) := / Ve|(c)-VC^.(c)dmO. (3.14c) 
Jn 



3.2.3 Finite dimensional case 



We assume now that the rate-of-shape-changes variable c has only a finite number 
of non-zero elements, say the N firsts {N > 1)- As explained in the Appendix, 
Subsection A.l, in this case we introduce the projector IIjv defined in (A. 2) and we 
identify Sn = I1n{S) with R^^. Upon this identification, the bilinear mappings 
N(c) and M'*(c) can be identified with actual matrices of sizes 3 x 2N and 2N x 2N 
respectively. Thus, we have: 



(c) := 



'In Ve[(c) • VC?(c) dmO Ve[(c) • Ve^(c) dm^' 

In ^^2ic) ■ Ve?(c) dmO Ve5(c) • Ve^(c) dm" 

In Va(c) • VC?(c) dmO • • • Va(c) • V^^(c) dm« 



while M''(c) reads: 



'(c) := 



/oVe?(c)-Ve?(c)dmO ... ^VC?(c).VC^(c)dmO 
/oV^?(c).VC?(c)dm^ ... ^VeJ(c).Ve^(c)d™0 

In Ve^(c) • Ve?(c) dmO . . . ^ Ve^V(c) • Ve^(c) dm^. 



XVe^(c).Ve?(c)dmO. . 

1/2 




/ove^(c)-ve^(c 









1/2 









1/iV- 




) dm!j 








1/A^ + l 



(3.15) 



(3.16) 



Note that in the finite dimensional case, the overall mass matrix M(c) can also be 
identified with an actual (3 + 2A^) x (3 + 2A^) block matrix defined by: 



1(c) := 



■M'-(c) 
N(c)^ 



N(c) ■ 
M^(c) 



(3.17) 
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3.2.4 Regularity of the Mass Matrices and Lagrangian function 

As a straightforward consequence of Theorem 3.1 we get: 

Theorem 3.2. e r{V,C2{'R.^ x R^)), N G r{V,C2i'R^ x S)) and M'* e 
V{V,C2{S X S)). It entails that M e V{V,C2{{^^ x 5) x (R^ x S))) and the 
Lagrangian function is smooth in all of its variables. 

Notice that the mass matrix M''(c) is the sum of M^(c) :— diag (m, m, /(c)) and 

"Lver-ve[dm/ ... /^vcr-vcsdm/- 

M^(c) : : 

inVeS-VCrdm/ ... /^Va-Vadm/ 

which is positive for all c e P. We deduce that, for all c e D: 

detM''(c) > detM5;(c) > m^TTy, (3.18) 

and therefore that M''(c) is invertible for all c E V. Further, according to the 
classical formula: 

(M'-(c))"' = ^^co(M'-(c))^, (3.19) 

^ ^ " detM'-(c) ^ ^ " ' ^ ' 

where co(M'"(c))^ stands for the transpose of the co-matrix, we deduce from esti- 
mate (3.18) and Theorem 3.2: 

Proposition 3.2. The application c E T) i-^ M'"(c)^-'^ E £(R'^,R'^) is well-defined 
and analytic with radius of convergence R > I. 

3.3 Explicit computation of the mass matrices 
3.3.1 Elementary potentials 

In this Subsection, our aim is to compute explicitly the elementary potentials defined 
in Subsection 3.1. According to (3.7), (3.9) together with expressions (3.10) and 
after some algebra, we get in polar coordinates: 

a„^[(c)(0, r) - cos(j^) + Jl^h (3-20a) 

a„^2^(c)(0, r) - cos(j-^) + J^o', sin(j0), (3.20b) 

a„^3(c)(6l, r)^Y^ kal cos{ke) + kal sin{k0), (3.20c) 



k>l 



and for k > I: 



dn^'kic){r,9) =Y,Jl^ijC0s{j9)+jfil^smUe), (3.20d) 

a„e^(c)(r,0) ^Y.^^h'^^'^W + J'^hMjO), (3.20e) 

where the sequences of real numbers {a^)k>i and (a^)fe>o are defined by: 

al — bj+iaj — flj+i^j, (3.20f) 

al^bk^i + '^bj+kaj - aj+kbj, (fc > 2), (3.20g) 
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and 



(3.20h) 
(3.20i) 



and the sequences {f^kj)j>i ^^^^ (mI- j)j>i = k> 1) by: 
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' {k/j + l)ak+j + {k/j ~ l)ak-j if 1 < j < fc - 1, 
{k/{k + 1) + l)a2fe+i - l/(fc +1) if j - fc + 1, 
[k/j + l)afc 



if = fc or J > fc + 2, 



{k/j + l)bk+, - (fc/j - l)bk-, iil<j<k-l, 



{k/j + l)bk+j 



iij>k, 



(3.21a) 
(3.21b) 



and 



[k/j + l)bk+j + {k/j - l)bk-j if 1 < j < - 1, 
(fc/j + l)bk+j if j > k, 

-{k/j + l)ak+j + {k/j - l)ak-j if 1 < j < fc - 1, 

-{k/{k + 1) + l)a2fc+i - l/(fc +1) if J - fc + 1, 

'{k/j + l)afc+j if J = fc or j > A: + 2. 



(3.21c) 
(3.21d) 



We deduce the following expressions for the elementary potentials defined in CI 
\ in polar coordinates: 



Ci{c){r,e 


) = E(^OjCos(j0) 


+ /^0jcos(i6»))r-^ 


(3.22a) 


C2{c){r,6 




^ 1^1^ cos{ je))r-\ 


(3.22b) 


C3{c){r,6 


) = ^(a]cos(j0) + 


a] cos{j9))r-\ 


(3.22c) 


Ck{c){r.6 




^ A,jC0s{j6))r'\ 


(3.22d) 


ek{c){r,0 




^vl^cos{j0))r-^. 


(3.22e) 



From now on, we will denote, for any fc > 0, /i.^, :— ((Mj; j)j>i: (/^fe j)i>i) P^i'^ of 
real sequences) and likewise := (('^fe j)j>i) ('^fc j)j>i) and a := ((a])j>i, (a^)j>i). 

3.3.2 Mass matrices 

The entries of the mass matrices defined in Subsection 3.2 can be now easily derived 
from the expressions (3.22) of the elementary potentials. Indeed, let us consider, 
for instance, the first element of the matrix M'"(c). Applying Green's formula, we 
get: 

/ VC[(c) • Ve[(c)dm/ = -pf j ^[(c)fl(c)da, 

Jo. JdD (JT 



and then: 



/ Vei:(c).Ve5:(c)dm/ = 7rp^.^j|Mo 
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For any two pair of sequences v :— ((uj)j>i, (w|)j>i) and ^ :— ((^j)j>i, (<r|)j>i) of 
real numbers, we introduce the notation: 

v-<;:=J2 H^l-^l + ^hl) and \v\^ ■.= v-v. 
fc>i 

Taking into account the expressions (3.12) and (2.15), it allows us to give the 
expression of the mass matrices in a convenient short form: 



'(c) = PoTT 



1 
1 





E 



\Ck\ 



2 ■ ^ fc+ 1 

fe>i 



(3.23) 



and likewise the elements (3.13) of the mass matrix N read, for all fc > 1: 



(pflTflQ-fl^ ifj=l 

[pf^^a■^l^. if j = 3, 

{PfTTfj.^ ■ fk a j = 1 

Pf^uo-vu ifj-2 

Pf-KOL ■ Vk if j = 3. 



(3.24) 



(3.25) 



At last, the expressions of the elements (3.14) of M''(c) reads, for all j, k > 1: 



^'^(c),a^a'=) = pfnfx^ ■ + (M'^(c), a^ b'=)(c) = pfirfx^ ■ uu. (3.26a) 



k + l 
k + l' 



(3.26b) 



In the finite dimensional case (i.e c has only a finite number of non-zero elements, 
the first, N > 1), treated in Paragraph 3.2.3, the expressions (3.15) and (3.16) 
turn out to be: 



N(c) = PfTT 



a ■ fii a. ■ Vi 



Mo • Mw Mo ■ 

- Mat 1^0 ■ t^N 
a - fij^ a - 



(3.27) 



and 



\c) := PfTT 



iMil Ml -1^1 
1^1 -Ml 



MjvMi Mat'^i- 



Ml Mat Ml ■ t^N 

1^1 t^N 1^1 t^N 



1/2 ... 
1/2 ... 






Trpo 



1/iV+l 
1/iV + l 



(3.28) 
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3.3.3 A special case, N = 2 

We specify N ~ 2 and we assume that both c and c have only two non-zero elements: 
ci = ai + ibi and C2 — a2 + ih2- The quantities arising in the expression of the 
matrix M''(c) are in this case: 

iMoP = (1 - ai)' + (&i)' + + 2(fe2)^ 

/^o • = -26i, 

/^o • a = 3(0261 - 0162) - 2aibia2 + 62[(ai)^ - (61)^], 

|i.ol'-(l + «i)' + (61)' + 2(02)' + 2(62)', 
I/O • a = 3(6162 + 0201) + 2016162 + a2[(ai)^ - (61)^], 
lap = [(ai)2 + (6i)2][(a2)2 + (62)^] + 2(61)^ + 2(01)^ + 3(62)^ + 3(a2)^ 



while the elements of the matrix N(c) read: 





Ml 


= 2(0201 + 6261) - 302, 


Mo 


t^l 


= 2(6201 - O261) 


-362, 


Mo 


At2 


= -Oi + (oi)2-(6i)2, 


Mo 




= —61 + 2oi6i, 






Ml 


= 2(0261 - 6201) - 362, 


1^0 


t^l 


= 2(6261 + 0201) 
= (6i)2 ~ (oi)2 - 


+ 302, 


1^0 


M2 


= 61 + 2oi6i, 


J^o 


1^2 


Oi, 


a. 


Ml 


= -6i-26i[(02)2 + (62)2], 

--62 + 62[(oi)2+(6i)2], 


cx 


t^l 


= ai + 2ai[(62)2- 

02 - 02[(6i)^ + 


(ain 


a. 


M2 


cx 


t^2 



Identity (2.5) together with (2.7) allow computation of the expression of the density 
p* in polar coordinates: 



pUx{c, r, 0)) = pa\l- (oi)2 - (61)^ - 4(0102 + 6i62)r cos(6i) 

+ 4(6102 - Oi62)r sm{e) - i[{b2)^ + {a2)^y 



Notice that this quantity is not required to compute the motion of the amoeba. 
Finally, to compute the internal forces of the swimming body, as it will be shown 
in Subsection 3.7.1, we need the expression of the elements of M'^(c). We give only 
the non-zero elements: 



2 ^4[(02)2 + (62)2] 

2(0261 -I- 0162), 



iMii 

Ml ■ 1^2 
fJ'2 '^1 =2(0162 + 6102), 
V1V2 ==2(6162-0201), 



Ml • M2 

IM2P 



2(0201 - 6261), 

(ai)2 + (6i)2 + 



3'_ 

= 4[(62)2 + (02)2] + 
W2\' =(6i)2 + (oi)2+i 



2' 



3.4 Equation of motion 

Following the method explained in [ , chap VI, pages 160-201], we introduce P 
and n, the translational and angular impulses, as well as L and A, the impulses 
relating to the deformations: 



Xc) 



and 



(N(c),c) 



In these identities, both left hand side terms can be identified with elements of R'^. 
We compute that for all p := (SjcD)-^ S R'^: 

ii(r*)-L -cj(s*)^' 




d d 


r 


d 


r 




dt 9q 








P = 
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We next easily obtain that: 



d^dL 
dt c)q 



P 



dL . _ d 
9q ^ dt 



P + L 




s* 




P + L 






n + A 




UJ 


+ 


n + A 








According to Theorem 3.2, the Lagrangian function is smooth with respect to all 
of its variables, allowing all of the derivatives to be computed. Invoking the least 
action principle, the Euler-Lagrange equation of motion is (see e.g. [21, Theorem 
7.3.3 page 187]): 



dt 9q 



9q 



0, Vp e r3 



In our case we get: 



-(n + A)-r*.(P + L)^ 



0, 
0. 



(3.29a) 
(3.29b) 



If at time t — Q: 



P + L 

n + A 



0, 



this relation remains true for all t > since the Cauchy-Lipschitz Theorem ensures 
the uniqueness of the solution of system (3.29). We obtain here the equation: 



-(M'-(c))-i(N(c),c). 



If we introduce, for alH > 0: 



T{t) :=rS+ / r*(s)ds, 
then we can rewrite (3.30) as a first order ODE: 

= -(M'-(c))-i(N(c),c), 



d_ 
di 



(t>0). 



(3.30) 



(3.31) 



(3.32a) 



This expression is the one obtained in [14] for articulated bodies. This expression is 
very convenient to study the motion of the shape-changing body since it gives the 
velocity with respect to the shape variable. Due to the change of variables (3.31), it 
has to be supplemented with a so-called reconstruction equation allowing to recover 
r knowing 6: 



r(i)=ro- / i?(0)(M'-(c))-i(N(c),c)ds. 



(3.32b) 



We can also easily give the equation of motion in terms of r and 6. To this purpose, 
we introduce the 3x3 block matrix: 



no) 



R{9) 
1 



and since r* = R{9)^r, we can rewrite (3.30) in the form: 



d 
dt 



-7^(6l)(M'■(c))-l(N(c),c) 



(t > 0). 



(3.33) 
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3.5 Mathematically allowable control function 



According to Definition 1, a control function t E [0, T] i-^ c(t) e 5 to be physically 
allowable has to satisfy the constraints c(t) G T), ||c(i)||7- — fi {jj, < 1) and identity 
(2.18b) for all t > 0. However, when c has only a finite number of non-zero elements 
(the N firsts, TV > 1), we observe that the expressions (3.23), (3.27) and (3.28) of 
the mass matrices make sense even if c ^ P and their entries are still polynomial 
functions in c. Likewise, the matrix M[''(c) is invertible for all c e Sn and the 
entries of M''(c)^^ are analytic functions in c, with infinite radii of convergence. 

However, when c g I?, the mappings x(c) and ^(c) may be no longer invertible 
and therefore the domains A* and J-* are ill-defined (they overlap themselves). 
The elementary potentials cannot be defined, either. But since we can consider 
expressions (3.23), (3.27) and (3.28) as defining abstract matrices (not relating any 
longer to our physical problem) for any c € Sn, we can also consider the ODE (3.33) 
in this case. It leads us to relax the constraint c e I? in the finite dimensional case 
and to introduce for all /i > and all integers > 0: 

£Ni^J') ■= {ce Sn ■■ ||c||r„ = A*}- 

We can next set: 

Definition 2 (Mathematically allowable control function). A continuous piecewise 
function t S [0,T] i— > c(t) G Sn (for some real positive T and positive integer 
N) is said to be mathematically allowable when: 

• There exists fi > such that c{t) G ^Nifj) for all t > 0. 

• Constraint (2.18b) is satisfied for allt g]0,T[ such thatc{t) is well-defined (as 
piecewise function, c is well-defined for all t g]0, T[ but a finite number). 







1 
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1 
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Figure 3: Examples of mathematically allowable shape-changes that are not physi- 
cally allowable (/i = 0.8, c ^V). The domain S* overlaps itself. The colors gives the 
value of the density inside the animal. Observe that we can have negative densities 
but the total mass of the animal is always positive and constant. 



3.6 Well-posedness 

From Proposition 3.2 and Theorem 3.2, we deduce: 

Proposition 3.3. For any fi in ]0, 1[, for any smooth physically allowable control 
function c : [0,T] — s- f*(/i) and for any initial condition (ro,0o) G R-^ R/27r, 
there exists one unique smooth solution to Equations (3.32) (or equivalently (3.33) j 
defined on [0, T]. 

We get the same result, replacing the physically allowable control by a mathe- 
matically allowable control function c : [0,T] £n{i^) for any integer N > 1. 

A mathematically allowable control function c is only assumed to be continuous 
and piecewise C^. It means that there exist to = < ti . . . < tn = T such that c be 

on each interval ]tk, tk+i[, fc = 0, . . . , n — 1. The solution given in Proposition 3.3 
is obtained by integrating the EDO (3.32) on each interval ]tk,tk+i[. It is also 
continuous, piecewise on ]0,T[. 
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3.7 Controlling with internal forces 



We have selected the shape of the amoeba to be the given quantity for controUing 
the motion of the animal. By defining the notion of physically allowable shape- 
changes (see Definition 1), we took care that these deformations result from the 
work of internal forces only. In this Subsection, we will first compute the expression 
of these internal forces in terms of the given shape-changes. Second, we will focus 
on the situation where the control variable c lives in the finite iV-dimensional vector 
space {N > 1) Sn (defined in the Appendix, Subsection A.l). We will prove that 
it is immaterial whether we select either the internal forces or the shape-changes 
as controls, the relation linking them being one-to-one. However, we will also show 
that the notion of allowable internal forces is much more difficult to define than the 
notion of allowable shape-changes. 



3.7.1 Expression of the internal forces 

First, the shape-changes being given, we are interested in computing the internal 
forces they result from. Let t S [0,T] i-^ c{t) G £'{^) (a* €]0, 1[) be any smooth 
allowable control in the sense of Definition 1 and compute, using the Euler-Lagrange 
equation (3.33) the corresponding induced rigid motion t £ [0, T] i-^ q(i) £ Q. The 
generalized forces, denoted by F in the sequel, are usually defined in Lagrangian 
Mechanics by: 

where L is the Lagrangian function defined in Subsection 2.5. This equality tells 
us that for all time, F is an element of S' (the dual space of S, defined in the 
Appendix, Subsection A.l). Let us then introduce the mass matrix K, defined for 
all physically allowable control c as an element of £2(5 x S) by: 

(K(c),6'i,c~2) := (M'^(c),ci,c~2)- (N(c),£i)^(M''(c))-i(N(c),£2), {ci,C2 £ S). 

According to Theorem 3.2 and Proposition 3.2, the mapping c £ V 1-^ K(c) £ 
£2(3 X S) is analytic. We next define the energy-like (or modified Lagrangian) 
amount: 

£:= i(K(c),c,c), {t>0). (3.35) 

One proves after some algebra and taking into account (3.33), that identity (3.34) 
can be rewritten as: 

„ d dC dC 

^-dtd-c-O^' ^'-''^ 
This expression allows us to compute straightforwardly the internal forces from the 
shape-changes without computing the induced motion q of the swimming animal. 
It can be slightly expanded. Since K(c) £ £2(5 x S), we deduce that, for all c £V: 

^(c) e £3(5x5x5), 
and the Frechet derivative of K in the direction c at the point c is given by: 



— (c),H,.,.; e £2(5x5). 

We deduce that the expanded form of (3.36) is: 

(K(c),c,.) + (r(c),c,c,.) = (F,.), (3.37) 
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where r(c) G £3(5 X 5 X 5) is a so-called ChristofFel symbol defined by: 



(r(c),Ci,C2,C3) 



/dK. , \ /dK. . \ /9K, , 

(^-T^(c),C2,Ci,C3 J + — (c),C2,C3,Ci J - — (c) , C3 , Ci , C2 



for all Ci,C2,C3 e S. In this form, F depends only on c and on its first and second 
derivatives with respect to time. The internal forces are the relevant quantities one 
has to consider when seeking optimal strokes. In [l.i] for instance, the cost function 
to be minimized over the time interval [0, T] is taken to be: 



J:= r\\F{t)\\l,dt. 
Jo 



3.7.2 Equivalence between controlling with shape-changes and with forces 

We assume now that the shape variable c lives in the finite dimensional vector space, 
Sn (for some integer iV > 1; see the Appendix, Subsection A.l). We use for the 
mass matrices the expressions (3.23), (3.27) and (3.28), allowing them to be defined 
for all c e Sn- Likewise, the mass matrix K(c) can be identified with an actual 
2A^ X 2N symmetric matrix: 

K(c) := M''(c) - (N(c))^(M''(c))-iN(c), 

and we claim: 

Lemma 3.2. There exists a constant vjq > such that, for all c, c G Sn: 

VK(c)e>^.A.||e|||„. (3.38) 

As usual, we refer to the Appendix, Subsection A.l for the notation. 

Proof. Let us set q :— (— M''(c)^^N(c)c, c)-'" and observe that: 

q'^M(c)q = c^K(c)c, 

where the matrix M(c) is defined in (3.17). We next easily get that, since M''(c) is 
positive: 

q^M(c)q> ^[M'^(c)-iN(c)c]'^M''(c)[M'^(c)-iN(c)c] + ^c'^M'^(c)c 

k=l 

where Ck — ak + ibk for all k G {1, . ■ ■ ,N}. The proof is completed after setting 
m = 2/ (pott N^). □ 

Remark 3.1. Observe that the conclusion of this lemma is no longer true in the 
general infinite dimensional case since vm as N oo. This explains why 
we are not able to prove the equivalence between controlling by shape-changes and 
internal forces in the general infinite dimensional case. 

The dual space of Sn can be identified with and equation (3.37) can be 
merely seen as an ODE in C^. We can now state our main equivalence result: 
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Theorem 3.3. For any given Lipschitz continuous function t E [0, T] t-^ F(t) G 
and for any Cauchy data (co,Co) € S'^r x 5jv, there exists a unique smooth maximal 
solution t e [0, T] i-^ c{t) G solving equation (3.37) and such that c(0) = Cq and 
c(0) = co. 

Remark 3.2. Here, an intricate (an up to now open) problem consists in deter- 
mining conditions for the generalized forces F ensuring that the shape-changes we 
obtain by integrating ODE (3.37) are allowable in the sense of Definition 2 (or even 
more complicated, in the sense of Definition 1). In other words, how to specify 
within Lagrangian formalism that the forces F are indeed internal to the animal? 

Proof. As explained in Subsection 3.5, all of the terms depending in c in the ODE 
are analytic. Further, Lemma 3.2 ensures that the matrix K(c) is always invertible. 
The Cauchy-Lipschitz Theorem applies and yields the existence and uniqueness of 
a maximal solution defined on some interval [0,T*) with T* < T. Remember that 
C = c-^K(c)c/2. One easily verifies that: 



d dL 
dt dc 



dC 



c = 



dC 

'dt' 



{te[o,T*)), 



meaning that the variation of energy C is equal to the power-like amount F • c. 
Integrating over [0,t] for any < t < T* and invoking inequality (3.38) we get: 



i^NWcfr^ < m < £(0) + / ||F(.s)||T„||c(s)||T„ds. 



(3.39) 



Setting then: 

T{t) := r|lF(,s)||Tjlc(.)|lT„ds, 
we obtain after some basic algebra that, for all < i < T*: 

r{t) ^ ||F(t)HT„ 
VT(<) + /:(o) - 

Integrating this inequality with respect to time, we get the estimate: 

_ o 

1 



T(i) < £(0) 



\F(s) 



ds 



2^A^ Jo 

Plugging this result into (3.39), we just have proved the Gronwall-type inequality: 



llTjv 



< 



-m 



1 

UN 



1 



I^N Jo 



t -|2 

l|F(s)||T„ds 



meaning that c remains bounded for all t G [0,r*). Classical behavior results for 
ODEs tell us that T* = T and the proof is completed. □ 



4 Controllability results 

This Section is dedicated to the study of control problems associated with Equations 
(3.32) or (3.33). 

4.1 Main Theorem of controllability 

We begin by giving our main controllability result. We make some comments and 
give the outline of the proof that will be next set out in the following Subsections. 
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4.1.1 Statement of the main theorem 



Theorem 4.1. For every fi^ in ]0,1[, for every e > and for every reference 
continuous curve (q^,c^^) : [0,T] — > Q x £'{ii^), there exists a real n in ]0,1[ and 
an analytic physically allowable curve c : [0,r] — > £*{^i) such that 

1. \\c{t) - {t)\\s < e for allt e]0,T[; 

2. The solution q of (3.33) starting from q^(0) satisfies ||q(t) — q^(i)||Q < £ for 
all t e]0,T[. 

Remark 4.1. Note that in the hypotheses of Theorem ^.l, we do not require that 
be physically allowable. More precisely, can violate constraint (2.12b). 

The conclusion of the Theorem may seem a httle bit surprising. It imphes that 
for given shape-changes t € [0,T] i-^ c^(i) g S, to which corresponds the motion 
t e [0,T] q^(t) e Q obtained by solving the Euler-Lagrange equation (3.33) (say 
for instance, moving forward), one can find shape-changes t S [0,T] i-^ c(t) G S 
arbitrarily close to t € [0, T] i-^- c^(i) e S (for the uniform norm on the set of maps 
from [0,T] to S), whose corresponding motion t e [0,T] i-^ q(t) G Q is arbitrarily 
close to any given trajectory (for instance, moving backward). This phenomenon 
will be termed Moonwalking. The strength of the Theorem can also be illustrated 
by selecting as reference function t ^ c^{t), shape-changes that do not result in 
locomotion. Such an example is given in the following paragraph. 

4.1.2 Flapping does not allow locomotion 

We establish in this paragraph a well known negative result (usually referred to as 
the scallop theorem, [29]). It states that it is impossible to achieve an arbitrarily 
large displacement of the amoeba by flapping. 

Let any physically allowable control function c : [0,T] ^'{lA be given (for 
some ^ e]0, 1[). 

Proposition 4.1. There exists a real number R > such that for any function 
f] : R+ — > [0,T] and for any initial condition qo E Q, the solution q^ :— (r^j^^)"^ : 
R+ — > Q of Equation (3.33) corresponding to the shape-changes t g R^. i— > c{P(t)) € 
£'{^), with initial condition qo remains in the ball of Q of center qo with radius R. 

Proof. Notice first that if c : [0,r] ^'(m) is physically allowable, then for any 
smooth function (3 : R+ — > [0, T], the control function c o /3 : i e R_|_ i-^ c{f3{t)) G 
£'{n) is also physically allowable. From the equation of motion (3.33), one gets 




(t) = -7^(0/3(/3(^)))M'■(c(/3(^))-l(N(c(/3(<))),c(/3(^)))/3'(<), 



which yields successively: 



q0(i) = q(0) - / n{ep{P{s)))M-{c{(3{s))-\N{c{(3{s))),c{p{s)))(3'{s)ds 



Jo 



= q(0)- / 7^(0^(T))M'^(c(r))-l(N(c(r)),c(r))dr. 




We next easily obtain the estimate: 




<r sup M'''(c(s))-i(N(c(s)),c(s)) . 



se[o,T] 



27 



The function t |M^^(c(s))N(c(s))c(s)| being continuous, it has a finite supre- 
mum A > on the compact set [0,T]. It entails that for every positive t, q(t) is 
contained in the ball centered in q(0) with radius R := TA. □ 

Remark 4.2. Proposition 4-1 does not apply any longer if the function [3 is no 
more assumed to he bounded. For instance, if c(0) — c(T) and (3 is the identity 
function, it may happen that t G R_|_ r^it) e is not bounded. See Section 5 
for examples of swimming with periodic deformations. 

4.1.3 Outline of the proof of Theorem 4.1 

The proof of the Theorem is somewhat intricate and will be done using finite di- 
mensional control techniques. In the following subsection, we state Theorem 4.2, a 
finite dimensional version of Theorem 4.1. Theorem 4.2 is proved in the particu- 
lar 2-dimensional case in Subsection 4.3 (using Lie brackets computed with Maple 
and Maxima, softwares allowing symbolic computations) and in the general case 
in Subsection 4.4 (the proof resting on the computations of Lie brackets done in 
Subsection 4.3). At last, in Subsection 4.5, we prove how the infinite dimensional 
problem of control can be suitably approximated by a finite dimensional problem 
for which Theorem 4.2 applies. This will conclude the proof of Theorem 4.1. 

4.2 Finite dimensional version of Theorem 4.1 

4.2.1 Statement of the Theorem 

Remember that we have denoted merely p the quotient po/pf > 0. 

Theorem 4.2. For every integer N > 2, for all but maybe a finite number of pairs 
{p, p) g]0, +cx)[x]0, +oo[, for every e > and for every reference continuous curve 
(q^,c^^) : [0,T] ^ Q X £n{p) (T > 0), there exists an analytic and mathematically 
allowable curve c : [0,T] — s- £n{p) such that 

1. \\c{t) - {t)\\s < e for allt e]0,T[; 

2. The solution q of (3.33) starting from q^(0) satisfies ||q(t) ~q^(i)||Q < £ for 
all t e]0,T[. 

Remark 4.3. As in Theorem 4-1, observe that the reference control function 
can violate constraint (2.18b). 

4.2.2 Restatement of the finite dimensional control problem 

In Theorem 4.2, it is required that the control functions t E [0,T] i-^ c{t) G Sn be 
allowable in the sense of Definition 2. It means in particular, according to (2.12b) 
and (2.18b), that (F(c),c) = (G(c),c) for all t e]0,T[ such that c{t) is well 
defined. In order to deal with these constraints, we consider Xjv (X^)i<j"<ti a 
set of n vectors fields in {n an integer, n > 1), where X-' := (X'J)fe>i and for all 
c e Sn, Xlic) xi{c) + iyl{c), x{{c), yl{c) e R, /c > 1 and = if fc > iV. 

Definition 3 (Allowable set of vector fields). A set of vector fields Xn :— (X-')i<j<„ 
defined in Sn is said to be allowable when: 

1. For all j G {1, . . . , n}, the field c G Sn ^''(c) G Sn is analytic. 

2. For allce Sn and all j G {1,.. .,n}, (i^(c), (c)) = (G(c), X^'(c)) ^ 0. 
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Consider now a set A :— (Aj)i<j<„ of piecewise constant functions from [0, T] C 
R into R and the EDO: 



c(i) = ^A,(t)X-'(cW), it>0). 



(4.1) 



The usefulness of Definition 3 arises with the following, easy to prove, property: 

Proposition 4.2. For any set of piecewise constant functions X = {Xj)i<j<n as 
above and any initial data Cq S Sn, there exists one unique solution t € [0,T] i— > 
c{t) G Sn to ED0{4:.1). Further, c is mathematically allowable (in the sense of 
Definition 2 with fi ||co||rjv/ 

Once the set has been chosen, the new control variables turn out to be 
the piecewise constant functions A = {Xj)i<j<n- We can give an example of such 
allowable set of vector fields when N = 2. Thus, let us define X2 (X.^)i<j<4 by: 



Xlic) = 

xKc) = 

xlic) = 
Xlic) - 



-2 [(6162/3 + aia2) + i (6102 - 0162/8)] , 
-2 [(0162 - 6ia2/3) + i (6162 + aia2/3)] , 
-(3/2) [(6162 + 0102/3) + I (0162/3 - 6102)] , 
-(3/2) [(6102/3 - 0102) + i (6162/3 -I- 0102)] 



Xlic) - |ci| 
Xlic) 
XHc) 
Xfic) 



l\Cl\ 

i\c2\ 



(4.2) 



where we recall that c := ick)k>i with Ck = Ok + ibk (ofe, 6^ e R) for all A: > 1. 
One can easily verify that X2 is indeed allowable in the sense of Definition 3. 
Going back to the general case, we can rewrite the EDO (3.33) in the form: 



d 
dt 



E"=lA,(^)7^(0)(M'-(c))-l(N(c),X^■(c))■ 



E;LiA,(t)x^-(c) 



(4.3) 



supplemented with initial conditions: (q(0),c(0)) — (qo,co) € Q x Sn- This then 
leads us to introduce the set of analytic vector fields 2)jv ■= (Y-')i<i<„, defined on 
Qx Sn by: 



Y^ (q,c) := 



-7^(6')(M'^(c))-l(N(c),XJ■(c))■ 
X^(c) 



(y^)l<j<n, 

V(q,c) eQxSN, 



(4.4) 



and to rewrite (4.3) as: 



d 
dt 



5^A,(t)Y^(q,c). 



(4.5) 



Notice that, according to the expressions of M''(c) and N(c), the vector fields 
depend not only on (^, c) but also on p :— po/pj, excluding all other quantities. 
Further, the dependence is analytic with respect to all of these variables (including 
p). Seen as a problem of control, equation (4.5) fits the general form of geometric 
control theory. So, before going further, let us now recall some important results. 



4.2.3 Tools of geometric control theory 

Let M be an analytic connected manifold, endowed with the Riemannian distance 
d-M and let A" be a set of analytic vector fields on M. 

Let (A'j)i<j<p (p G N, p > 1) be a finite sequence of elements of X. For any 
given finite sequence (i/)i<;<p of positive real numbers, we define T :— Ym=i ™^ 
the application 

'i'(t,A),<,<, :Mx [0,T]->M, 
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by: 

*(t,A)i<,<p('7,0) = q 

This definition can easily be extended to the case of non-complete vector fields, but 
the domain of $ is then restricted to the product of a certain (possibly empty) open 
set of M by the interval [0, T]. 

Definition 4.3. We say that the trajectories of X can track a given continuous 
curve 7 : [0,r] M {T £ R_|_) if for every e > 0, there exists a finite sequence 
S = {tj, Xj)i<j<p of elements of R x A" such that (i) 7(0) x [0,r] belongs to the 
domain of 4>5, (ii) $5(7(0), T) = c(T) and (iii) dM($s(7(0), t), 7(t)) < e for all 

te [0,T]. 

The following proposition is a pretty classical consequence of the Orbit Theorem 
(recalled in Section B). We give a proof for the sake of completeness (see [10] for 
further discussions). 

Proposition 4.3. If X is a symmetric cone such that liieqX = TqM for any q in 

M, then the trajectories of X can track any given continuous reference curve on M . 

Proof. Let e > be a positive number and let 7 : [0, T] — > M be a continuous curve 
on M. Since [0,T] is compact, 7 is uniformly continuous on [0,T]. Hence, there 
exists some 77 > such that for any t, t' in [0, T], \t — t'\ < 77 implies dM{j{t),'^{t')) < 
e/3. Without loss of generality, we may assume that T — Nrj with G N. For 
any integer n in [0,N — 1], define M„ as the open e/2 neighborhood of j{nri), 
and consider Xn the restriction of X to Af„. The set M„ is a connected analytic 
manifold, which contains both ^{nrf) and 7((n + 1)»7). For every q in Mm the Lie 
algebra of Xn at q is equal to TqMn- From Proposition B.4, we deduce the existence 
of a finite sequence S'„ — (i", X")i<;<p^^ such that <^s,^ iliil ^ ^)''i)jV) = 7(^^); and 
$S^(7((Z — l)r]),t) belongs to M/ for every t in [0, 77], / in {1 . . .p„}. 

Define now S as the concatenation of So, S2, ■ ■ ■ , Sjy^i. The set S* is a finite 
sequence of elements of the product of R and X. By construction, {7(0)} x [0,T] 
belongs to the domain of $5 and $5(7(0), T) — 'y{T). For every t in [0,r], there 
exists an integer n not greater than N — I such that nrj < t < (n + I)?]. Hence 
$5(7(0), t) = $5^(7(n?7,t - nifj)) belongs to A/„, that is 

dM($5(7(0), i), 7W) = dM($s„7('^'7, i ~ n7]),-/{t)) 

< dM{^s„l{nr],t ~ nT]),-f{nri)) + dM{7{riri),-f{t)) 
e e 

and the proof is completed. □ 

The main idea in the proof of Theorem 4.2 will be to apply Proposition 4.3 on 
the analytic manifold Q x Sn and with analytic set of vector fields 2} at defined in 
(4.4). 

4.2.4 Remark on the optimal control problem 

Geometric control theory gives for free the existence of optimal solutions to our 
control problem. In particular, Filipov Theorem (see [1, Chapter 10]) ensures that: 

Theorem 4.4. Let f : Qx iSjv x R" —^'R.bea continuous function, IC be a compact 
subset o/R", (qo, Cq) and (qi, Ci) be two points ofQxSN- If there exists a trajectory 
of (4.5) steering (qo,co) to (qi,ci), associated to a measurable bounded control A 
taking value in IC, then there exists a measurable bounded control : [0,T] — > R" 
taking also value in IC such that: 
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1. The corresponding solution (q^,c^) of (4.5) steers (qo,Co) to (qi,Ci); 

2. The triplet (q^,c^^,A^) realizes the infimum of the cost 

''^ f{q{t),c{t),X{t))dt, 

among all measurable bounded controls A taking value in K, and steering the 
system (4.5) from (qo,co) to (qi,ci). 

Notice that if there exists a trajectory linking (qo,Co) to (qi,ci), then Cq and 
Ci both have to belong to the same set £Ar(/i) for some /i > i.e., they have to 
satisfy ||co||r„ = l|ci||rjv = M- 



4.3 The case = 2 

In this section, we focus on the case N — 2 and we consider the set of analytic 
vector fields X2 defined on £2(m) (for all ^ > 0) by (4.2) . As explained in the 
Appendix, Subsection A.l we identify £2{fJ') with E2{fi) the 3-dimensional analytic 
submanifold of R*. Furthermore, the dynamic induced by the vector fields X2 on 
S2{fJ') is the same as the one induced on E2{fi) by the vector fields X2 ■= {X^)i<j<i 
defined for any c := (ai, 61, a2, 62)"^ S by: 



X\c) 



-2(6i&2/3 + aia2)' 
-2(6ia2 - 0162/3) 



if + bi 




05 



bl 







-3(6162 + aia2/3)/2 
-3(ai62/3- 6102)72 



X'ic) 



X\c) 



-2(0162 - 6102/3) 
-2(6162 + 0102/3) 


al + bl 







-3(-oia2 + 6i02/3)/2 
-3(6162/3 + oia2)/2 



It is easy to verify that for every c S R^, c ^ 0, the linear space spanned by X2 
has cardinal two. Since the expressions of the mass matrices M''(c) and N(c) are 
given in Subsections 3.3.2 and 3.3.3 we can compute explicitly the expression of the 
vector fields in 2)2 := (Y^)i<j<4 defined in (4.4) in Q x 52. Like X2, 2)2 can be 
identified with the set of vector fields 3^2 '■= (^■')i<j<4 defined in Q x R"*. 

4.3.1 Computation of the Lie algebra 

If we do not require the control function to be analytic in Theorem 4.2, according 
to Proposition 4.3, it suffices to check whether Lie(3^2) has dimension 6 everywhere 
on Q X i?2(M) to prove the theorem when N — 2. Since the expressions of the fields 
(1 < J < 4) are somewhat intricate, we first concentrate on the set Lie(A2). 

Proposition 4.4. For any fi > 0, the family X2 is completely nonholonomic on 
E2{n), that is, for any c in R*, c 7^ 0, Liec(A2) ~ TcE2{fi) where ^ := ||c||t2- 

Proof. Recall that for every non-zero c e R^, the four vectors X^{c), X'^{c), X^{c) 
and X*{c) span a 2-dimensional subspace of the 3-dimensional linear space TcE2{n), 
{fi := IIcIIts)- a direct computation (see Proposition B.2) gives 



[X\X']ic) = -{al 



bl) 



-bl' 
oi 
-362 

302 



and [X^,X'^]{c) 



{4 



bl) 



-bl 

Ol 

-362 
3a2 
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For any c = (01,61,02,62) such that of + bi ^ 0, the vectors -'^^(c) and X'^{c) 
are clearly linearly independent. Proceed by contradiction and assume that X^{c), 
X'^{c) and [X^,X^](c) are not linearly independent. Then, there exists ai and a2 
(two real numbers) such that [X^,X^](c) = aiX^{c) + a2X'^{c). Since a\-\-h\ ^ 0, 
one has ai — —462 and a2 — 4a2. Projecting the equality [X^,X^](c) = aiX^{c) + 
a2X'^{c) on the first coordinate, one gets — 6i(ai + h\) = 26i(a2 + 62), that is 



02 = ^2 = or 61 = 0. If 02 = 62 = 0, then 



6^ = which is in contradiction 



with 



h{ ^ 0, hence 61 = 0. From the projection on the second coordinate. 



one gets (of + 6f)ai = — 2ai(a2 + 63), from which one deduces ai = 0, which is 
incompatible with the hypothesis a\ + b\ —Q. We have then proved that for any c 
such that a\ + h\^ 0, -^^(c), ^^(c) and [X^,X^](c) are linearly independent. 

The same argument shows that, if 0^+63 ^ 0, then X^(c), X^c) and [X^ , X^c) 
are also linearly independent. 

For every non-zero c € R**, dimLiec(<%'2) > 3 = dimTc-B2(Ai)- Therefore, 
Liec(A2) = TcE2{fJ.) and the proof is completed. □ 

Since the family X2 is completely non holonomic on every submanifold E2{^), 
the attainable set at any positive time of the control system 



c{t) ^ \j{t)X\c{t)), c(0) = Co e 



with piecewise constant controls {^j)i<j<A is equal to E2{^i) (/i := ||co||t2)- 

We next consider the control system defined as the projection of the system 
(4.3) on R/27r x ^2. Indeed, denoting by {v)^ the third component of any vector v 
of R'^, this system reads 



d 
dt 



c 

c)(0) = 



■-E,tiA,(t)(M'-(c)-i(N(c),X^(c))^^ 
E-=iA,(t)X^(c) 

9o,co) eR/2^xf2(M)■ 



We define X2 (X.^)i<j<4 the set of analytic vector fields (!<.? < 4) in 
R/27r X 52 by 



(m'-(c)-1(N(c),XJ'(c)))_ 
X^-(c) 



, (1 < i < 4). 



According to the expressions of the matrices 1^(0) and N(c), the vector fields X^ 
(1 < j < 4) depend on c and also on p po/pf, excluding all other quantities and 
the dependence is analytic. As previously, the set X2 can be identified with the set 
of vector fields A2 := (X-')i<j<4 defined in R/27r x R'*. 

Proposition 4.5. For all but maybe a finite number of pairs (/i, p) g]0, +oo[x]0, +oo[ 
and for any (6, c) G R/27r x E2{p), the Lie algebra Lie(-g c)('^2) is equal to the whole 
tangent space T(^q ,.)(R,/2'k x E2{p)). 

This proof involves symbolic computations performed with both software pro- 
grams: Maple and Maxima. The relevant computations can be downloaded on 
the web page http://www.iecn.u-iiancy.fr/~munnier/page_amoeba/control_ 
index.html. 

Proof. The proof comprises three steps: 

Step 1: Let be /x > and p > Q and choose any point c* e E2{p). Since the 
vector fields composing X2 do not depend on 9, we deduce that the Lie algebra 
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Lie(e.c*)('^2) has the same dimension d, (0 < < 4) for any 9 G R/27r. According 
to Proposition 4.4, the orbit of X2 through c* is equal to £'2(/^) and hence the orbit 
of X2 through any point {Oi,c) E R/27r x i?2(M) contains at least one point {62, c*) 
for some 6*2 G R/27r. As a consequence of the Orbit Theorem (Theorem B.2), 
Lie(0j c)('^2) is also of dimension d and hence, the dimension of Lie^g (.)('^2) is con- 
stant on R/27r x E2{fJ,)- 

We deduce that for any fixed pair (fi, p) e]0, +oo[x]0, +oo[, it is enough to prove 
that there exists at least one point (6'*,c*) G R/27r x £2(1.1) for which Lie(g. c.)(A2) 
has dimension 4 to prove that Lie(0 c)('^2) = r(0 c-)(R/27r x E2{^)) for all {6,c) G 
R/27r X E2{p). 

Step 2: Explicit computations of the iterated Lie brackets of the fields are 
straightforward, yet quite intricate. In order to simplify the resulting expressions, 
we define the fields M = det(M'~)XJ and — (^i)i<j<4. Since p'^'^ det{yr) 
does not vanish, the linear spaces Lie(g and Lie(e have the same di- 

mension (see Proposition B.3). Further, like the fields , the fields Xi depend 
only (and analytically) on c and p :— p^/pf. Despite this simplification, the ex- 
plicit results of the computation are definitely still too long to be printed. Since 
the fields Xi do not depend on 0, these expressions also only depend on c and p. 
We chose a specific c having the form c_|_ := (a;,0,x,0) and c_ := (— x, 0,— x, 0) 
(a; G R). When x describes R+, p := ||c+||t2 = ||c_||t2 = VSa; describes also R+. 
We use Maxima and Maple to obtain the huge expressions of Xl{6, c+), X^{9, c+), 
[Xl,Xj]i0,j,+ )j,ml [Xl, [Xl,X^]]{9,c+) and Xl{0,c_), X^{9,c^), [XlX^]{e,c^) 
and [Xl, [Xl, X'^]]{9,C-) in terms of x and p. With each family of 4 column vec- 
tors, we define a 5 x 4 matrix and extract a 4 x 4 submatrix. Computing next 
the determinants of these square matrices, we obtain two polynomials P'^{p,x) and 
P~{p,x) of degree 20 in the variable p and 65 in the variables x. We first consider 
them as polynomial in the variable x with coefficients in R(p), the field of rational 
fractions in the variable p with coefficients in R. Seeking their gcd, we obtain a 
monomial. If we denote by Zp the finite set consisting of all of the zeros of either 
the numerator or the denominator of all of the rational fractions of R(p) arising in 
the Euclidean algorithm leading to the expression of the gcd, we deduce that for all 
p ^ Zp, P'^{p,x) and P^{p,x) (seen as polynomials in the variable a:) have only 
as common root and next that for any a; ^ 0, at least one out of the two families of 
vectors fields spans a 4-dimensional vector space. This proves the Proposition for 
all of the pairs {p,p) g]0, -|-oo[x (]0, -|-oo[\Zp). 



P 




Figure 4: In red, the possibly fead pairs {p,p) of ]0, -|-oo[x]0, -|-oo[ 

Step 3: Let us now speciahze x = 1 {p = a/3) and see both P~^{p, 1) and P~{p, 1) 
as polynomials in p with coefficients in R. Their gcd (also computed with Maple 
and Maxima) is equal to 1, meaning that they have no common root. This yields 
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the conclusion of the Proposition for all of the pairs (\/3, p), p > 0. Let us next pick 
some p^ £ Zp and consider {P^{p^,x) — P^{p^,x))^ + P~^{p^,x)^ as a polynomial 
in X. If it vanishes in an infinite number of points, then it would be identically 
zero, which would contradict the result obtained for x = 1. We deduce that for any 
p g Zp, the polynomials in x, P~^{p,x) and P~{p,x) have at most a finite number 
of commune zeros and the proof is completed. □ 

Remark 4.4. The Authors conjecture that the result of Proposition 4-5 actually 
holds true for all of the pairs [p, p) e]0, +oo[x]0, +oo[. 

We denote *B the finite (maybe empty) set of all the pairs [p, p) e]0, +cxd[x]0, +oo[ 
for which the conclusion of Proposition 4.5 does not hold true and we claim: 

Proposition 4.6. For all but maybe a finite number of pairs {p, p) e]0, +oo[x]0, +oo[ 
and for any (q, c) G Q x E2{p), the Lie algebra Lie(q c)(3^2) is equal to the whole 
tangent space T(q c)(Q x E2{p)). 

Proof. The proof is very similar to that of Proposition 4.5, so we only give the 
outline: 

Step 1: Since the expressions of the vector fields in 3^2 do not depend on r, the 
same arguments as those of step 1 in the proof of Proposition 4.5, allow us to prove 
that for all (p, p) E (]0, +oo[x]0, +oo[) the dimension of Liejq c)(3^2) is constant 
on Q X E2{p)- We deduce that for any fixed pair {p, p) G (]0, +oo[x]0, +oo[) \ 03, 
it is enough to prove that there exists at least one point (q*,c*) G Q x E{p) for 
which Lie(q.^c*)(3^2) has dimension 6 to prove that Lie(q^c)(3^2) — r'(q,c)(Q x E2{p)) 
for all (q,c) S Qx E2{p). 

Step 2 and Step 3: We define Y,^ := p-^^ detWY^ for all 1 < j < 4 and denote 

y* '■— (i^."')i<j<4. Using the same notation as in the proof of Proposition 4.5, 
we consider the points (qo,c+) and (qo,c_) with qo := (0,0,0)-^ G Q. Using 
Maxima and Maple, we next jprove that for all but maybe a finite number of pairs 
(/i, p) G (]0, +oo[x]0, +oo[) \ 05, at least one out of the two families: 

{y;(qo,c+),r.2(qo,c+),[y.i,y.2](qo,c+), 

[Y,\ [y.\y.2]](qo,c+), [y,\yMcio,c+), [[y,\y,% [y,\yM^„,c+)} 

and 

{ y.l (qo , c_ ) , ( qo , c_ ) , [1.1 , y.2] (qo , c_ ) , 

[Y,\ [r;,r.2]](qo,c_), [y:,r.2]](qo,c_), [[y,\y,% [y.\y.2]](qo,c_)}, 

spans a 6-dimensional vector space. □ 

We define 05 as the (possibly empty) subset of ]0, +oo[x]0, +oo[ (containing 05) 
that is excluded in the statement of Proposition 4.6. 

4.4 Proof of Theorem 4.2 

We perform now the proof of Theorem 4.2 in the general case. 
4.4.1 A distribution on £n{p) 

Fix N an integer {N > 2) and let us define Xat a set of n :— 2N{N — 1) vector fields 
on Sn- The vector fields composing Xn are denoted X^'^^-'^i^i) := {xj!"'''^^''^^)k>i 
and -X^koM.^) {xi''°'''''^^)k>i with 1 < ko ^ h < N. For any couple (fco,fci) 
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of distinct integers less than N, all of the components of the fields X'^'^^''^^''^^ and 
X('=o,fci,2) g^j,g constant equal to zero, but (maybe) the four following entries: 

Xi1'"'-'^^k,{k, + l){al + bl), 

-^(ko,ku2) ^ (^-klakobki - kiakobki + k^ak^bko + kaak^bko) 

+ i {-kjbk„bk-, - kibk„bk^ - klak„aki - fcoafc^afci) , 
xj^"''^^'^=^koik, + l)ial + bl). 

One can verify that Xat is allowable in the sense of Definition 3. 

Remark 4.5. The fields Xt^'^^i), X(i'2^2)^ x'^2a,i) andX.^'^'^^'^^ extend to Sn (N > 
2) the fields X^, X^, X^ and defined in S2 by relations (4.2). 

More generally, for all integer N' such that 2 < N' < N , we have Xpf G ^nIsj^i ■ 

As in Paragraph 4.2.2, we introduce ^n, the set composed of the vector fields 
Yfco:feij defined for all (q, c) in Q x En{ij) by: 



-7^(6')(M'■(c))-l(N(c),X('^0''=l■J')(c))' 
X(feo,feij)(c) 



for j — 1,2 and 1 < feg 7^ fci < and we rewrite (4.5): 

n 

- E A(fc„,,,,,)(t)Y('=o.fe...)(q^e), (4.6) 



dt 



i = 1.2 



where X^ko.kij) '■ [^tT] — > R are piecewise constant functions. 

Remark 4.6. Like in Remark 4-5, for all integers N' such that I < N' < N , we 
have C ^nIs^, ■ 

As already mentioned in Subsection 4.3, via the identification between Sn and 
Sn := R^^, we can identify Xn with Xn, whose elements are denoted X'^'^"''^!'^^ 
and j = 1,2, 1 < fco 7^ A;i < A^. The system (4.6) turns out to be a 

finite dimensional control system on the analytic 2N + 2-dimensional submanifold 
Q X i?Ar(^) in R^^+'^ (as usually, :— ||co||tni where (qo,co) G Q x Sjy is the 
initial condition of System (4.6)). 



4.4.2 Computation of the Lie algebras 

The following Proposition is a generalization of Proposition 4.4 to any dimension 
A > 2. 

Proposition 4.7. For all N > 2 and for any /i > 0, the family Xn is completely 
nonholonomic on En{ij), that is, for any c in R^^, c ^ 0, Liec(A'jv) = TcEn{ij) 
where fi := ||c||t„ . 

Proof. Let c be as in the Proposition. Since c is not zero, there exists fcg such that 
^ko +^feo ^ ^' "^^^ vector space span{A(*"'''^'-'^(c) : k =^ ko, j = 1,2} has dimension 
2Af — 2. To prove that the Lie algebra Liec(Ajv) has dimension 2A^ — 1, it is enough 
to check that (^jv(c), A('=0''=1'2)]) ^ g j^^. i^^^ ^^j^g mapping 

Fn being defined in the Appendix, Subsection A.l). A simple computation shows 
that {FNic),[X('^oMA)^x('^oM.2^) ^ _2^p^^^+alJk^{k2 + l)iblk, + alk, + blko + 
al^ko) < for any choice of fci 7^ fco, and the proof is then completed. □ 
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As in Paragraph 4.3.1, we next define on R/27r x Sn the set of vector fields Xjv 
which can be identified with defined on R/27r x Sn- According to our usual 
notation, we denote X^'^"-''^-^'' (j = 1, 2, 1 < /cq 7^ fci < A^) the n elements of Xn- 
We can restate Proposition 4.5 in the general A^-dimensional case: 

Proposition 4.8. For all integer N > 2, for all of the pairs {fi, p) e (]0, +oo[x]0, +oo[)^ 
*8 and for any {0,c) e R/27r x E^di), the Lie algebra Lie(e c)('^Af) is equal to the 
whole tangent space T(e c)(R-/27r x En{^j))- 

Proof. When N = 2, the result is given by Proposition 4.5 so let us assume that 
A^ > 3. As in the proof of Proposition 4.5, we first establish that for all N > 
3 and all {fi, p) g]0, +oo[x]0, +oo[, the dimension of Lie(e^c)('^A') is constant on 
R/27r X ENifJ.)- We next define <% := {A^^^'^i'^), 1 < j < 2, 3 < /ci < A^} C AV. 
Observe that for all A^ > 3 and for all {p-,p) e]0, +oo[x]0, +oo[, f2(Ai) C Swit^) 
and hence E2{p) can be seen as an immersed submanifold of En{p). We can 
easily verify that for all {0,c) £ R/27r x E2{p) {E2{p) seen as a subset of Ejq(ix)), 
dim(span(g j.^(A'jJ^)) — 2{N — 2) and Licjg (.)(A'2) n span^g — {0}, whence we 

deduce that: 

Lie(e,c)('^2) + spaujg = Lie(e_c)('^2) © span(e,c) (-^A) C Lie(g^e)('^7v)- 

But in Proposition 4.5, we have proved that for all of the pairs {p, p) E (]0, +oo[x]0, +00 
S and for any (0,c) g R/27r x E2{p), dim(Lie(g^c)('^2)) = 4. We deduce that 
dim(Lie(6).c)('^7v)) > 2A^, which is the dimension of R/27r x En{p), and the proof 
is completed. □ 

We can also generalize Proposition 4.6 as follows: 

Proposition 4.9. For all N > 2, for all pairs {p, p) € (]0, +oo[x]0, +oo[) \ 53 
and for any (q, c) G Q x En(p), the Lie algebra Lie(q c)(3^Ar) is equal to the whole 
tangent space T'(q.c)(Q x Epf{p)). 

Proof. The proof is similar to the proof of Proposition 4.8, reusing the results of 
Proposition 4.6. □ 



4.4.3 Regularization 

According to Proposition 4.3 together with Proposition 4.9, we have proved so far, 
that for any integer N > 2, for any pair {p,p) G (]0, +oo[x]0, +oo[) \ *B, for every 
e > and for every reference continuous curve (q^^,c^^) : [0,r] Q x £n{p), there 
exists 2A^(A^— 1) piecewise constant functions (XkQ.ki.j) i<fco^'=i<«. • [0) T] R, such 

that the solution (q,c) of (4.6) starting from (qt(0), ct(0)) reaches (q^(T), ct(r)) 
at time T and remains e-close to the reference curve (q^, c^) for all time between 
and T. 

The proof of Theorem 4.2 follows, because the analytic real functions are dense 
for the L^{[0,T]) norm in the set of measurable bounded functions. One can there- 
fore approximate the piecewise constant control functions {Xj)ko,ki-j (1 ^ ^0 7^ ^ 
A", j = 1, 2) by a suitable family of analytic functions. 



4.5 Proof of Theorem 4.1 

4.5.1 Finite dimensional approximation 

Proposition 4.10. Let e > 0, p := pn/pf > and p^ G]0,1[ be given and : 
[0,T] — > £'{p^) be a continuous curve (not necessarily physically allowable). Then, 
there exist p g]0,1[ such that {p,p) G (]0, -|-cx)[x]0, -|-cx)[) \ *B, an integer N > 2 and 
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a continuous curve c 

[o,r]. 



[0,T] £n{^J-) such that \\c(t) — c'(t)||5 < s for every t in 



Proof. For every integer N > 2 and every e' > define O^r := {t G [0, T] : 
\\c^t) — IlNC^{t)\\s < £'} (the projector Hjy from S onto Sjy is defined in the 
Appendix, Subsection A.l). Because is continuous, the set Oat is open in [0,T] 
for all A'^ > 2 and since for any t e [0,r], UNC^t) ct(t) as A^ ^ oo, we 
deduce that [0,T] C Un>iQn- The interval [0,T] being compact and the sequence 
{Qn)n>i non-decreasing, [0,T] C 9^ for some A^ (depending on e'). However, we 
cannot yet choose 11 ^yc^ as a good finite dimensional approximation of because 
||n7vc^'(t)||7- is certainly not constant. Since ||c||r < ||c||5 for all c g 5, we have: 



~ e' < \\IlNc\t)\\T < Vte[0,T]. 



(4.7) 



We deduce that \\UNC^t)\\T ^ for all t e [0,T] if e' is chosen small enough. For 
any /i in ]0, 1[ and for for any t in [0, T] we have next the estimates: 



njv(ct(t)) 
'l|n^(ct(i))lh 



ct(t) 



< 



Ml 



nAr(ct(t)) 

nA.(ct(i))|| 



r 



+ ||n^(cT(<))-cT(t)| 



< 



<c 



niv(ct(<))||r 

^N{cHt))\\T 



1 

- 1 



where C := maxt(=[oT] ||njv(c^(t))||5. According to (4.7), we deduce that, for all 
tG[0,T]: 



M] 



njv(ct(^)) 
n^(ct(i))||. 



ct(i) 



< C max • 



^ — /i ' I I ^ - _|_ 



t - 



-1} 



■ £ . 



(4.8) 



For any £ > 0, one can always choose s' small enough and /i satisfying the require- 
ments of the Proposition, such that the right hand side of (4.8) be smaller that e. 
We conclude the proof by choosing c := /in7v(c^)/||nAr(c^)||7-. □ 



4.5.2 Conclusion of the proof 

Let fi^ e]0, 1[, £ > and a reference continuous curve (q''',c''') : [0,T] — > Q x 
be given. Apply next Proposition 4.10 with e/2 to obtain an integer N > 2, fi e]0, 1[ 
and a continuous curve : [0,r] i-^ £jv(a') such that ||c"f(t) — c^(f)|j5 <e/2. Finally, 
apply Theorem 4.2 with £/2 and reference curve (q^,c^) : [0,r] i-^ Qx£fq{^), to get 
the conclusion of Theorem 4.1. Observe that since we are able to find c : [0,r] i-^ 
£N{y) C £{^J) such that \\c\t) — c(t)||5 < e for all £ > and since c^{t) G V for all 
t e [0,T], we can always assume that c is valued in £*{fi). 



5 Numerical results 

By integrating equations (3.32a) and (3.32b), we can easily compute the trajectory 
of the swimming animal. Indeed, all of the mass matrices 1^(0) and N(c) arising 
in the ODE have been made explicit in Subsection 3.3. 

This Section is accompanied by a web page containing further animations and 
numerical experiments and is located at: http : //www . iecn . u-nancy . f r/~munnier/ 
page_ainoeba/control_irLdex.html. We will always choose pf = 1 (the density of 
the fluid) and the density of the amoeba will be next computed accordingly based 
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on formulae (2.14) (neutrally buoyant case). The values of /i depend on the dimen- 
sion N of Sn (the Banach space of the control variable c) . They are chosen in such 
a way that £^{^i) = ^NifJ-) (i-e. the ball of center and radius /i of T/v be included 
in V). All of the animations and figures have been realized with MATLAB. 

5.1 Swimming using Lie brackets 

Although recourse to Lie brackets can be useful to derive theoretical controllability 
results, they yield in general a quite inefficient swimming strategy. This can be 
illustrated by the following example. We consider here the configuration described 
in Subsection 4.3 i.e. we specialize N = 2, ^ = 1/2, Tq = (0,0)-^, = and 
Co = (1/2,0)"^. The shape-changes and the trajectory are given by integrating 
the EDO (4.3). We approximate the displacement induced by the Lie bracket 
by integrating the EDO with first Ai = 1 and Xj = for j = 2,3,4 over 
a small time interval (of length 0.1) then we set X2 ~ I and Xj = for j — 1,3,4 
over a time interval of the same length, next Ai — —1, Xj = for j — 2,3,4 
and finally A2 = — 1 and Xj — for j — 1,2,3. We repeat this process fifty 
times to obtain the trajectory of the center of mass (i.e. the parameterized line 
r{t) := {ri{t) , r2{t)) with t G [0,20]) of the amoeba pictured in Figure 5 while in 
Figure 6 we display the shape variables ci and C2 with respect to time over the time 
interval [0, 20]. A movie related to this simulation is given on the web page http: 
//www . iecn . u-nemcy . f r/~munnier/page_amoeba/ control_index . html, showing 
even more clearly how inefficient this swimming strategy is. 

0.06 



0.05- 
0.04- 



0.03- 
0.02- 
0.01 - 

0- 

-0.02 0.02 0.04 0.06 

Figure 5: Trajectory of the center of mass of the amoeba when the shape-changes 
are obtained by approximated Lie brackets. The bold (red) line can be considered 
as the net displacement of the swimming animal. 
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Figure 6: Shape variables ci{t) = ai{t) + ibi{t) and C2{t) — 02 (i) + 162 (^) with 
respect to time. At time t = 20, the hnes correspond respectively from top to 
bottom to ai{t) (red), 62 (^) (black), 02 (t) (blue) and bi{t) (green). 

5.2 Examples of more efficient swimming strategies 

Throughout this subsection, we specify TV = 6 (i.e. the shape variable reads c = 
(ci,...,C6)^ e with Cfe := Ofc + ibk for all 1 < fc < 6). The first difficulty we are 
faced with in seeking swimming strategies is that the control variable 1 1-^ c{t) has 
to be allowable in the sense of Definition 1. To reflect this constraint, we introduce 
the new functions 1 1-^ aj{t) (1 < j < 5) and 1 1-^ hk{t) (1 < < 3) that we choose 
to be our new control functions and we next define: 



Ri{t) :— ^cos{ai{t)) 

R2{t) 



sin(Q!i(i)) cos(a2(i)), 

v2 



i?3(i) := 



: sin(Q!i(i)) sin(a2(t)) cos(Q!3(t)), 



Ri{t) :— ^ sin(Q!i(<)) sin(Q!2(^)) sin(Q;3(t)) cos(Q!4(i)), 

Rz{t) :— sin(Q!i(i)) sin(Q:2(t)) sin(a3(i)) sin(Q!4(t)) cos(a5(i)), 
v5 

i?6(^) '■— sin(Q!i(i)) sin(Q!2(t)) sin(a3(i)) sin(Q!4(t)) sin(Q;5(t)), 
V6 



and 



h{t) 



3 7o 
1 

5 Jo 



hi{s)Rl{s)ds, 



hi{s)Rl{s)ds, 



1 



/i2(s)i?4(s)ds, 04{t):^ 7 / h2{s)Rl{s)ds, 



h3{s)Rl{s)ds, 96{ty.= 



4 7o 



1 



6 



h3is)Rl{s)ds. 



If we set now: 

ak{t) := Rk{t) cos{0k{t)) and bk{t) := Rkit) sm{dk{t)), 
then it can be readily verfied that the function 1 1— > c(t) is indeed allowable. 
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5.2.1 Example of a straight-forward motion 

We set ^l = 0.5, ai{t) = t, aj{t) = (j = 2,3,4,5) and hk{t) = (fc = 1,2,3) for all 
t > to obtain a net straight- forward motion for the amoeba consisting of periodic 
strokes (of 2tt time period). With these data, the functions 1 1-^ Ck{t) {k — 3, 4, 5, 6) 
are constant equal to zero for all t > 0. Screenshots of the amoeba over a stroke 
are given in Figure 7 while in Figure 8 are drawn the x-coordinate of the center of 
mass of the animal and the x-coordinate of its velocity with respect to time. 
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Figure 7: Screenshots of the motion of the amoeba over a stroke. The colours give 
the value of the internal density. The animal is neutrally buoyant, so at rest its 
density is 1 (the density of the fluid) . 
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Figure 8: x-coordinate of the center of mass of the amoeba (dashed line) and x- 
coordinate of its velocity (solid line). 



Following the method described in Subsection 3.7.1, the shape-changes being given, 
we can compute the expression of the internal forces. Like the shape-changes, the 
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internal forces are also 27r-periodic as it can be seen in Figure 9. 
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(a) Controls ai{t) (solid line) and 02 (t) (dashed line) with 
respect to time over a stroke. Both controls f)i(t) and 62 (t) 
are equal to zero. 
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1 2 3 4 5 6 

(b) Coordinates Fi(t) (solid line), Fi{t) (dashed line) and 
F4{t) (dashed-dot line) of the generalized internal force over 
a stroke. The coordinate F2{t) is equal to 0. 



Figure 9: Shape changes and internal forces. 



5.2.2 Example of circular motion 

We set now again ^ = 0.5, ai{t) = t, aj{t) = {j — 2,3,4,5) and hk{t) ~ 
{k — 2,3) for all i > but we specify the control variable hi to be a non-zero 
constant. With these settings, we observe that the amoeba swims along a circular 
path. 




-3-2-101234 

Figure 10: Successive positions and shapes of the amoeba in its course when hi = 1. 
The animal follows a circular trajectory completed over a time interval of length 
approximately 247r. 
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Again, with these data, only ci and C2 are non-zero functions. The radii of the 
circles change along with the values of the constant hi as illustrated in Figure 11. 
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Figure 11: Trajectories of the center of mass of the amoeba. For hi — —1.5 in red 
over the time interval [0,187r], for /ii = 1 in green over [0,247r] and for hi = 2 in 
blue over [0, 127r]. 

When hi — I, the graphs of the controls are given in Figure 12 and the graphs 
of the internal forces in Figure 13. 




-0.51 ' ' ' ' ' ^ 
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Figure 12: Values of the controls ai{t), bi{t), 02 (i) and 62 (^) over the time interval 
[0,247r] for the amoeba following the circular trajectory with hi ^ 1. 
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Figure 13: Components Fi{t), F2{t), F3{t) and F4{t) of the generalized force over 
the time interval [0, 247r] for the amoeba following the circular trajectory with hi = 
1. 

5.2.3 Motion planning 

Based on the two preceding examples, we observe that the amoeba can follow any 
smooth trajectory with only the two first control variables ci and C2 being non- 
zero. Indeed, the function ai governs the frequencies of the strokes (and hence the 
velocity of the animal) and hi can be seen as the steering function; it allows the 
amoeba to turn left or right. On the web page, http://www.iecii.u-iiancy.fr/ 
~inunnier/page_ainoeba/control_iiidex.html, such examples of motion planning 
are given. 
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5.2.4 Further examples 



The preceding example in fact illustrates a more general feature of locomotion: 
motion planning is actually possible with any pair of consecutive controls variables 
{ck,Ck+i) {k > 1). For example, if we are willing to use only the pair (03,04) as 
controls, then we can set ai{t) — a2{t) = tt/2, Q4(<) = a5{t) = 0, hi{t) — /i3(i) — 
for all i > and the frequencies of the strokes are driven by while ft,2 turns out 
to be the steering variable. Such a strategy is illustrated in Figure 14. 



1 



0- 



-1 - 



-2 



-3 



-4- 




-5- - 

-3-2-101234 

Figure 14: Successive positions and shapes of the amoeba for ai{t) = ct-iit) = 7r/2, 
a4{t) = Q5(t) = 0, hi{t) = hsit) = for alH > and asi^t) = t, h2{t) = 1.2. 

At last, we can use the pair (05, cg). We set then ak{t) — tt/2 {k = 1,2,3,4), 
hi{t) — h2{t) = for all < > and the variable 0:5 controls the frequency of the 
strokes while /13 becomes the steering variable. These settings are those used in 
Figure 15. 




-2.5- 



-8 -6 -4-2 2 

Figure 15: Successive positions and shapes of the amoeba for ak{t) = 7r/2 {k = 
1,2,3,4), hi{t) = h2{t) = 0, a5(t) = t and hsit) = 1 for t e [0,1007r], h^it) = -2 
for t e]1007r, 2007r], h^lt) = for i e]2007r, 3007r] and h^^t) = 4 for < e]3007r, 4007r] 

Once more, we refer to the web page http : //www . iecn .u-nancy . f r/~munnier/ 
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page_Eimoeba/control_index.litml for the animations and further examples. 

5.3 Moonwalking 

Our main result of controllability, Theorem 4.1, states that the amoeba is not only 
able to follow approximately any given trajectory but also that this task can be 
achieved while undergoing (also approximately) any prescribed shape-changes. Let 
us illustrate this surprising result with the following example: In Figure 16 are 
displayed screenshots of the amoeba swimming to the left in the first row and 
toward the opposite direction in the second row although the shape-changes seem 
to be similar in both cases. 
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Figure 16: On the first row, the animal swims forward while it swims backward on 
the second row. The shape-changes seem to be similar in both cases. 



Obviously, the shape-changes are actually different but the difference cannot be 
observed at a (time and space) macro-scale. The controls are set to be in both 
cases: ai{t) = t, a2{t) = tt/6, a3(t) = 7i'/12, a4(t) = 7r/12 and hk{t) — 
{k = 1,2,3) for all t € [0,47r] but a^it) = in the first case while a^it) = -lOOOOt 
in the second one. It means that in the second case and with respect to the 
first case, there is superimposed shape-changes with very high frequency and very 
low amplitude making the animal swimming backward. We refer to the web page 
http : //www. iecn.u-ncincy . f r/~munnier/page_amoeba/ corLtrol_index .html for 
relating animations and further explanations. 



A Additional material relating to the modeling 

A.l Banach spaces of complex series 
A. 1.1 Infinite dimensional spaces 

Remember that we denote the complex sequences c := {ck)k>i with = Ofe + ibk, 
P = —1 and Ofc, 6fe S R for all fc > 1. All the sequences we consider in this article 
hve in the spaces: 



S := {(cfe)fe>i e C'^ : ^fc(|afc| + \bk\) < +oo}, 

k>l 

T :- {(cfc),>i e : 5] fc(|afep + \b,\') < +oo}. 



k>l 



which are Banach spaces once endowed with their natural norms. We denote B the 
unitary open ball of S and 



V (c := (cfc)fc>i e 5 : sup V(fc + l)cfc+iz'' < l|. 

zedD I > 



k>0 
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According to the inequality ||c||7- < ||c||5 for all c E S, we deduce that S C T. 
If we define = (a;^)fe>i and b-' = (6^)fe>i G S {j > 1), the complex sequences 
such that a^. = 61 (the Kronecker symbol) and bl. = iSj., then {a^, h^,j > 1} is a 
Schauder basis of T and S. Using the Hilbert structure of T to express the duality 
product, the dual of S can be identified with: 

S' := {(cfe)fc>i e : max{|afe|, |5fc|, fc > 1} < +c»}. 

We introduce next: 

c ^ ^^(c) c ^ G(c), 

where, for all c := (cfe)fegN, = afc + ife^, we set: 
(F(c),c) := ^ ^-^(flfefefc - fefcOfc) and (0(0), c) := ^ ^(0^0^ + 6^5^). (A.l) 

fe>l k>l 

A. 1.2 Finite dimensional spaces 

The projector n^r is defined for any integer > 1 and any complex sequence 
c = (cfe)fc>i by: 

. ^^ fcfc if fc < A^, , . , 

njv c fe = < " A.2 

^ ^ " [0 iffc>Ar. ^ ' 

We denote Sn ■— ^n{S) and T/v ■— njv(T). Throughout the paper, we will identify 
Sn and T)v with, respectively, Sn and T/v that are nothing but (or R^^), but 
endowed respectively with the norms: 



N 



fc=i 



N 



Y.H\a,\'+\b,\') 



.k=l 



1/2 



for all c := (cfe)i<fc<Ar £ . We recall the standard inequalities, for all integer 
A > 1: 

||c||t„<||c||s^<v/7V(7vTT)||c||t„, (ceC^). (A.3) 

Remember that for all integers A^ > 1 and all /i > 0, we have defined Sm{ij) '■= {c G 
Sn ■ l|c||r = /i}- As for the space Sn and Sn, we have identified fAr(^) with 

N 

EN{y) := {(01,61,02,62, • • .,aN,bN) € R^^ : ^ k{al + bl) = ,u^}, 

k=l 

which is diffeomorphic to the 2N — 1 dimensional Euclidian sphere. For any c, c € 
C^, the definition (A.l) of F and G leads us to introduce: 

(F7v(c),c) := ^ — Y(afe6fc-6A;afc) and {Gn{c),c) ■.= ^k(akak + bkbk)- (A. 4) 

k=l k=l 

Remark A.l. According to (A.3), when /i < 1/ ^yN{N + 1), En{pi) and En{^j) 
are both included in the unitary ball of their ambient spaces. In particular, for such 
values of N and fj, and for any control function t G [0,r] c(<) G £N{f) there 
is no difference between being Physically allowable in the sense of Definition 1 and 
being Mathematically allowable in the sense of Definition 2. 
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A. 2 Neumann boundary value problem 

Note that this section is self-contained and independent, including the notation. We 
recall some results about the well-posedness of a Neumann boundary value problem 
in an exterior domain. 

Let be the exterior of a compact in R^, assume that H. is Lipschitz continuous, 
connected and consider the following general problem: 

— Am = in ri, dnU = g on dU, (A. 5) 

where g is a given function in L^(dfl) satisfying the so-called compatibility condi- 
tion: 

/ <?da = 0, (A.6) 

and n is the unitary normal to dfl directed toward the exterior of f2. The compat- 
ibility condition leads us to introduce Lj^{d^l) := {g S L^{d^l) s.t. (A.6) holds}. 
We denote by 'D'{^1) the space of distributions, define the weight function p{x) := 
[■y/l -t- jxp log(2 -I- and introduce the quotient weighted Sobolev space: 

Hl^{n) := {V- e V'{n) : piP £ L^{n), d^^^/j e L^n), Vz = 1,2}/R. 

The quotient means that two functions differing only by an additive constant are 
equal in this space. It is a Hilbert space once endowed with the scalar product and 
associated norm: 

(w,«)Hi,(n) ■■= J^Vu-Vvdx, \\u\\Hj^(n) — ")h^(o)- 

The variational formulation of (A. 5) is: 

/ Wu-Wvdx= / gvda, Vw S (A.7) 
Jn Jan 

and Lax-Milgram Theorem ensures that this problem admits one unique solution 

Theorem A.l. The linear operator: 

■.g£L%idn)^ueHlf{n), (A.8) 
where u is the solution of (A.7), is well-defined and continuous (i.e., A^^ g 

ciL%idn),Hl,m). 

A. 3 Regularity of the potential function 

According to (3.9) and (3.10), we can state that for all z G dD and c E V: 

a„4"(c)(z) := bl^+{bl„c)iz), a,4(c)(z) := b^^ + {bl„ c) (z) , 
where we have defined for all z G dD: 

6^,o(z) := -5R(z'=+i), (6^,^, c)(z):^^ + _ fca^, 

bl,{z) := -3(z'=+i), (6^,1, c)(z) -.^Y. + (^) - ^^fc- 

j>i-k ^ 

It is clear that for any fixed c G all of the series converge normally in z on dD 
and hence their sums define continuous functions on dD. However, we will rather 
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consider them as functions of Lj^{dD). For all c € I?, c = (ofe + ibk)k>i € S and 
z G dD, we have the following decomposition of the boundary data of (9„^''(c), c): 

k>l 

In this form and based on simple estimates, we can easily prove that dnS,'^ € 
P{'D,C{S,Lj^{dD))). We invoke next the continuity of the linear operator A]^^ 
defined in (A.8), to conclude that ^'^ € V{V,C{S,Hl;{n))). 

B A brief survey of the Orbit Theorem 

hi this Appendix, we aim to recall the statement of the Orbit Theorem used in Para- 
graph 4.2.3 to prove some trackability properties. The material presented below is 
now considered as a classical part of geometric control theory. It has been intro- 
duced in the beginning of the 20th century independantly by Rashevsky (1938) and 
Chow (1939) following the ideas of Caratheodory (1909). The results were unified 
and generalized by Jurdjevic and Sussmann in the 70's. In the following, we have 
chosen a very simplified presentation (restrained to the symmetric analytic case) of 
the exposition of [']. Many other very good textbooks present this material with 
detailed proofs and comments, see for instance [9], [11], [26] or the research papers 
[34] and [12]. 

Throughout this section, M is a real analytic manifold, and Q a set of analytic 
vector fields on M . We do not assume in general that the fields from Q are complete. 

B.l Attainable sets 

Let / be an element of Q and q* be an element of M . The Cauchy problem 

4 = /(?), 9(0) -<Z*, (B.l) 

admits a solution defined on the open interval I{f,q*) containing 0. For any real 
t in /(/, q*) we denote the value of the solution of (B.l) at time t by e^^ {q*). We 
denote by /(/, g*)+ = /(/, q*) n ]0, +oo[ the positive elements of /(/, q*). 

For any element q^ in AI and any positive real number T, we define the attainable 
set at time T oi Q from qo by the set Aqg{T) of all points of M that can be attained 
with Q using piece wise constants controls in time T 

Aq„{T) = {e*-^^ o e'^-i^"-! o • • • o e'^f'iqo) : p e N, /, e 

U e /(/„ e*-i-^-i o . . . o e*i/i(qo))+, ti + • • • + = t}, 

the times ti and the fields f i being chosen in such a way that every written quantity 
exists. We define also the orbit of G trough qo by the set Oqg of all points of M 
that can be attained with G using piecewise constant controls, at any positive or 
negative time 

AoiT) = {e*--^- o e*"-!^"-! o . . . o 6*1-^1 (go) : peN, f.eG, 

£/(/,, e*'-i-^'-io... e*i-^^(<Zo))}- 

Of course, if C/ is a cone, that is if Xf G G for any positive A as soon as / belongs to 
G, the set Aqg (T) does not depend on the positive T but only on go- G is assumed 
to be symmetric, that is if — / belongs to G as soon as / belongs to G, then the 
orbit of G trough a point go is the union of all attainable sets at positive time of G 
from go. 
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B.2 Lie algebra of vector fields 

If /i and /2 are two vector fields on M and g is a point of M, the Lie bracket 
[fii f2\{q) of /i and /2 at a point q is the derivative at < = of the curve t ^ 7(v^) 
where 7 is defined by ^{t) := e~*^^e~*f'-e*^^e^^'-{q) for t small enough. The Lie 
bracket of /i and /2 at a point q is an element of the tangent space TqM of M 
at the point q. The Lie bracket is bilinear and skew-symmetric in fi and /2, and 
measures the non-commutativity of the fields /i and /2 (see [ , Prop 2.6]). 

Proposition B.l. For any fi, /2 in Q, we have the equivalence: 

for all times ti and ^2 (tf ony) for which the expressions written in the left hand 
side of the above equivalence make sense. 

Lie brackets of vectors fields are easy to compute with the following formulas 
(see [1, Prop 1.3] and Exercise 2.2]). 

Proposition B.2. For any f\, /2 in Q. for any q in M , 

[A, /.](.) = f.M.)-f /.(.). 

Further, we have the useful property: 

Proposition B.3. Let fi and /2 be two smooth vector fields on M , and let a, h : 

M ^ be two smooth functions. Then 

[aX, bY] - ab[X, Y]+ ( ^x] Y-( ^y] X. 



^dq J \dq 

From the Lie brackets, we can define the Lie algebra: 

Definition B.l. The Lie algebra of Q is the linear span of all Lie brackets, of any 
length, of the elements of G 

Lie g = span{[/i, [. . . [f^-u fk] ■ ■ ■]], k G N, /, e G}, 

which is a subset of all the vector fields on M. 

We denote by LieqG ■— {g{q), g & Lie G} the evaluation LiCqG of the Lie algebra 
generated by at a point q oi M. 

B.3 The Orbit Theorem 

The Orbit Theorem describes the differential structure of the orbit trough a point 
(see for instance [I , Th 5.1] for a proof). 

Theorem B.2 (Orbit Theorem). For any q and qo in M : 

1. O{qo) is a connected immersed submanifold of M . 

2. Lfqe Oiqo), then TqOiqo) = Uc^G . 

Remark B.l. The conclusion (1) of the Orbit Theorem holds true even if M and 
G are only assumed to be smooth (and not analytic). The conclusion (2) is false in 
general when G is only assumed to be smooth. 

The Orbit Theorem has many consequences, among them the following useful 
properties (see [1, Th 5.2] for a proof and further discussion). 

Theorem B.3 (Rashevsky-Chow). Lf LieqG = TqM for every q in M, then the 
orbit of G through q is equal to M. 

Proposition B.4. LfG is a symmetric cone such that Lie^^ ~ TqM for every q in 
M , then the attainable set at any positive time of any point of M is equal to M . 
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